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ABSTRACT 

'Hie  axisymmetric  nonsimilar  compressible  laminar  boundary- layer 
equations  including  approximate  transverse  curvature  terms  were  modified 
to  treat  inviscid  external  vorticity,  slip  and  temperature  jump  as 
first-order  quantities,  and  this  is  referred  to  as  the  first-order 
treatment.  The  effects  of  boundary-layer  displacement  were  also 
treated.  Primary  interest  was  in  predicting  experimentally  measurable 
quantities  over  the  entire  body  length  of  nonanalytic  shapes,  and  the 
analysis  was  not  confined  to  the  nose  or  stagnation  region.  A  review, 
is  included  of  second-order  boundary-layer  theory  and  of  recent 
developments  in  the  numerical  solution  of  second-order  boundary- layer 
effects.  Comparisons  of  predicted  displacement-induced  pressure 
distributions,  heat-transfer  distributions  and  zero-lift  drag  were  made 
with  results  from  second-order  boundary-layer  theory  and  experimental 
data  for  a  spherically  blunted  9-degrec  half-angle  cone  at  free-stream 
Mach  numbers  of  9  and  18.  At  moderate  to  high  Reynolds  numbers,  the 
comparisons  showed  good  agreement  between  first-  and  second-order 
predictions  and  experimental  results  for  drag  and  heat-transfer  dis¬ 
tributions;  however,  poor  agreement  was  found  between  predicted  and 
experimental  displacement-induced  pressure  distributions.  At  low 
Reynolds  numbers,  both  first-  and  second -order  treatments  substan¬ 
tially  overprcdicted  the  zero-lift  drag.  The  range  of  applicability 
of  the  theories  was  established  for  the  conditions  treated  by 
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inspection  of  the  numerical  results  and  by  comparison  of  the  numerical 
results  with  experimental  zero-lift  drag  data.  Separate  and  coupled 
first-  and  second-order  boundary- layer  effects  were  considered  for 
experimentally  measurable  quantities.  The  first-order  treatment  showed, 
for  example,  that  slip  and  temperature  jump  produced  a  small  effect 
when  considered  as  a  separate  effect;  however,  this  effect  became 
important  and  under  certain  conditions  dominant  when  coupled  with 
other  effects  at  low  Reynolds  numbers.  Limitations  and  difficulties 
in  both  first-  and  second-order  boundary-layer  treatments  are  dis¬ 
cussed. 
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NOMENCLATURE 


Unless  otherwise  noted,  all  lengths  are  nondimensionalized  by  the  nose 
radius . 


sh 


cDr 


'disp 


'LC 


\STJ 


Base  area 

Parameter  defined  by  Equation  40f 

Nose  radius  in  second-order  boundary- layer  theory 

Slip  constant, 

Parameter  in  displacement  thickness  with  external  vorticity 
(see  Appendix  C) 

Nose  bluntness  parameter,  rn/ry, 

Nose-to-shock  radius  ratio,  rn/rsh 
Slip  constant,  3/4 

Density  viscosity  product  ratio,  py/p  u 

c  6 

Friction-drag  coefficient  referenced  to  base  area 
Pressure-drag  coefficient  referenced  to  base  area 
Displacement  term  in  continuity  equation  defined  by 
Equation  38a 

2 

Skin-friction  coefficient f  2 xw/p  U 

9  W  OO  CO 

Longitudinal  curvature  term  in  continuity  equation  defined 
by  Equation  38b 

Slip  and  temperature  term  in  continuity  equation  defined 
by  Equation  38a 
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CpyC  Transverse  curvature  term  in  continuity  equation  defined  by 
Equation  38c 

Cvort  Vorticity  term  in  continuity  equation  defined  by  Equation  38a 
Coefficient  in  linear  viscosity  law,  (Pw/O  (T«/Tw) 

Cj  Slip  constant,  (15/8)  (ir/2) 

Cp  Constant  pressure  specific  heat 

Edisp  Displacement  term  in  energy  equation  defined  by  Equation  40d 

Elc  Longitudinal  curvature  term  in  energy  equation  defined  by 

Equation  40b 

ESTJ  Slip  and  temperature  jump  term  in  energy  equation  defined 

by  Equation  40e 

ETVC  Transverse  curvature  term  in  energy  equation  defined  by 

Equation  40c 

Evort  Vorticity  term  in  energy  equation  defined  by  Equation  40a 
e  Error  test  in  first-order  momentum  and  energy  equation 

solutions 

Fs  Parameter  defined  by  Equation  86c 

f  Dimensionless  stream  function 

f1  Velocity  ratio,  u/ue 

g  Stagnation  enthalpy  ratio,  H/He 

H  Stagnation  enthalpy 

h  Static  enthalpy 

hre£  Reference  specific  enthalpy,  2.119  x  10®  ft2/sec2 

j  Index  0  or  1  for  two-dimensional  or  axisymmetric  flow, 

respectively 
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K 

k 

L 

t 

M 

Mdisp 

MLC 

mstj 

^rvc 

**vort 

N 

!. 

n 

P 

Pr 

P 


Parameter  defined  by  Equation  85b;  bound  on  $'(n  )  in  first- 

00 

order  treatment  (see  Appendix  A) 

The  factor  1  +  «n 

Characteristic  length  (dimensional) 

Axial  length  of  body  (dimensional) 

Mach  number 

Displacement  term  in  tangential  momentum  equation  defined 
by  Equation  39d 

Longitudinal  curvature  term  in  tangential  momentum  equation 
defined  by  Equation  39b 

Slip  and  temperature  jump  term  in  tangential  momentum 
equation  defined  by  Equation  39e 

Transverse  curvature  term  in  tangential  momentum  equation 
defined  by  Equation  39c 

Vorticity  term  in  tangential  momentum  equation  defined 
by  Equation  39a 

Inner  expansion  variable  in  second-order  boundary- layer 
theory,  n/e;j  parameter  in  first-order  treatment  defined 
below  Equation  78 

"Outer  expansion  variable  in  second-order  bound ary- layer 
theory,  n*/a* 

2 

Inviscid  outer  flow  pressure,  p*/p  *U* 

CO  00 

Prandtl  number 

Pressure 

Free-stream  normal  shock  pitot  pressure 
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Q  Total  velocity  vector 

R  Inviscid  outer  flow  density;  radius  parameter;  gas  constant 

Re  Reynolds  number,  p^U^*  rn*/y* 

Re«  i  Reynolds  number, 

Re*  Reynolds  number,  p*u*  rn/u* 

r  Radius 

rn  Nose  radius  (dimensional) 

S  Entropy 

Sj  Entropy  derivative  in  basic  inviscid  flow,  dS^/df 

St  Stanton  number,  q*/[p*U*c* (TQ*-TW*) ] 

s  Surface  distance  measured  from  stagnation  point 

T  Inviscid  outer  flow  temperature,  T*/(U*2/c*) 

“  P 

T  First-order  transverse  curvature  parameter 

t  Inner  flow  temperature 

U  Tangential  component  inviscid  outer  flow  velocity 

Free-stream  velocity 
u  Tangential  velocity 

V  Normal  component  inviscid  outer  flow  velocity 

v  Normal  velocity 

v^  Hypersonic  viscous  interaction  parameter,  M^CC^/Re^  ^ ) ly^2 

x  Surface  distance  from  stagnation  point 

y  Distance  normal  to  surface 

z  Axial  distance  from  stagnation  point 

a  Angle  between  tangent  to  surface  and  axis  of  symmetry 

0  Pressure  gradient  parameter,  (x/ue)  due/dx 
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Y 
6 

6* 

6* 

D 

6* 

v 

e 

t 

n 

n. 

6 

c 

K 

X 

y 

p 

o 

T 

* 

♦ 

Y 
* 
fi 

w 


Ratio  of  specific  heats,  1.4 
Boundary- layer  thickness 

Boundary- layer  displacement  thickness,  6**/a* 

Davis*  displacement  thickness,  6*/e 
Displacement  thickness  with  external  vorticity 
Van  Dyke's  expansion  parameter  defined  by  Equation  8 
Inviscid  vorticity 
Transformation  variable 

Edge  of  boundary  layer  in  transformed  plane 

Cone  half  angle 

Longitudinal  curvature 

Second  viscosity  coefficient 

Dynamic  viscosity  coefficient 

Density 

Prandtl  number  in  second-order  boundary- layer  theory 
Shear  stress 

Dissipation  function  defined  by  Equation  5 
Transformation  variable,  f-n 
Outer  stream  function 

Transformation  variable,  g-1;  inner  stream  function 
Vorticity  index  defined  by  Equation  79 
Van  Dyke's  inviscid  vorticity 
Exponent  in  power  viscosity  law 
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SUBSCRIPTS 

o  Stagnation  conditions 

1  First-order  quantity 

2  Second-order  quantity 

b  Base 

disp  Displacement  effect 

e  At  the  edge  of  the  boundary  layer 

eff  At  the  effective  body  surface 

i  lnviscid 

LC  Longitudinal-curvature  effect 

N  Derivative  with  respect  to  inner  variable  N 

n  Derivative  with  respect  to  outer  variable  n 

sh  Shock  wave 

STJ  Slip-and-temperature-jump  effect 

TVC  Transverse-curvature  effect 

vort  Vorticity  effect 

w  Wall 

x  Derivative  with  respect  to  x-coordinate 

y  Derivative  with  respect  to  y-coordinate 

*  Sonic  conditions 

“  At  free-stream  conditions 

SUPERSCRIPTS 

1  Derivative  with  respect  to  independent  variable 

*  Dimensional  quantity  in  second-order  theory 
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(k)  Iteration  index  in  first-order  displacement  treatment 

In  the  figures,  open  symbols  and/or  curves  are  used  to  denote  numerical 
results  and  filled  symbols  denote  experimental  data. 
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CHAPTER  I 

INTRODUCTION  AND  DEFINITION  OF  TERMS  USED 


Before  the  introductory  discussion  of  the  problem  is  presented, 
a  few  terms  will  be  defined  which  may  not  have  a  clear  definition  in 
the  literature.  After  the  terms  are  defined,  the  background  to  the 
problem  will  be  given  and  the  methods  used  and  scope  of  the  research 
reported  in  the  dissertation  will  be  described. 

I.  DEFINITION  OF  TERMS  USED 

Entropy  layer.  This  is  an  inviscid  layer  of  gas  near  the  wall 
which  has  passed  through  the  curved  portion  of  the  bow  shock  wave.  The 
entropy  in  this  layer  decreases  from  a  maximum  along  the  wall  to  a 
minimum  along  a  streamline  passing  through  the  inflection  point  in  the 
bow  shock  wave  (if  one  exists)  and  increases  to  the  constant  value 
bohind  a  straight  (conical)  shock  wave;  This  concept  will  be  made 
clearer  later  in  this  chapter  by  consideration  of  the  inviscid  shock 
layer  over  a  spherically  blunted  cone  at  hypersonic  flight  conditions. 

First-order  boundary-layer  effects.  Prandtl’s  classical 
boundary-layer  equations  for  conservation  of  mass,  momentum  and  energy 
for  two-dimensional  and  axisymmetric  flow  aro  the  basic  first-order 
equations.  If  terms  are  added  to  this  system  of  equations  to  approxi- 
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mate  other  effects  or  if  the  boundary  conditions  for  the  equations  are 
modified,  the  resulting  system  of  equations  is  called  first  order, 
and  the  numerical  solution  method  is  called  a  first-order  treatment  of 
higher-order  boundary- layer  effects.  It  is  possible  only  a  posteriori 
to  determine  the  magnitude  of  any  specific  effect  by  comparison  of 
solutions  with  and  without  the  effect  included.  Moreover,  in  general, 
simple  scaling  laws  do  not  exist  to  scale  wall-measurable  quantities 
such  as  skin-friction  coefficient  or  heat-transfer  rate;  therefore,  a 
series  of  calculations  is  necessary  over  the  range  of  the  variables  of 
interest. 

Second-order  boundary- layer  effects.  If  the  flow  variables  such 
as  velocity,  temperature,  density,  pressure,  etc.  are  expanded  in  terms 
of  an  appropriate  small  parameter  and  the  resulting  expansions  are 
substituted  into  the  Navier-Stokes  equations,  it  is  possible  to  obtain 
two  systems  of  equations  describing  the  flow.  The  lowest-order  system 
is  Prandtl's  equations,  which  are  called  the  classical  first-order 
boundary- layer  equations.  These  equations  are  nonlinear  partial  differ¬ 
ential  equations.  A  second  system  of  linear  partial  differential 
equations  is  obtained  where  various  terms  can  be  identified  with 
second-order  boundary- layer  effects  of  physical  significance.  The 
second-order  effects  are  not  all  independent  and  uncoupled,  but  with 
one  exception  simple  scaling  laws  can  be  obtained  to  scale  the  results. 

Longitudinal  curvature  effect.  When  the  surface  of  a  body  is 
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curved  rather  than  flat  in  the  direction  of  the  free-stream  velocity- 
vector,  centrifugal  forces  produce  a  pressure  change  normal  to  the 
surface.  This  curvature  will  influence  both  the  inviscid  and  viscous 
flow  fields;  however,  the  effects  on  the  inviscid  flow  field  are  usually 
larger,  persist  further  downstream  along  the  body,  and  are  more  easily 
treated  numerically  than  are  the  effects  on  the  viscous  boundary  layer. 

Transverse  curvature  effect.  A  second  curvature  effect  arises 
when  axisymmetric  bodies  of  revolution  are  considered.  The  boundary 
layer  on  a  body  of  revolution  grows  less  rapidly  than  the  change  in  wall 
radius  with  axial  distance  far  downstream  along  the  body  because  of  the 
transverse  curvature  effect.  Moreover,  the  radial  distance  of  points 
in  the  boundary  layer  from  the  axis  of  symmetry  may  not  be  assumed  equal 
to  the  local  wall  radius  when  the  ratio  of  the  thickness  of  the  boundary 
layer  to  the  local  wall  radius  cannot  be  assumed  negligible. 

Displacement  effect.  In  classical  boundary- layer  theory,  the 
properties  at  the  edge  of  the  boundary  layer  are  obtained  from  the 
inviscid  flow  field  solution  over  the  geometric  body  without  a  boundary 
layer.  The  boundary  layer  retards  the  flow  and  displaces  the  stream¬ 
lines  near  the  surface.  This  effect  on  the  inviscid  outer  flow  appears 
as  the  flow  over  the  geometric  body  thickened  by  the  boundary-layer 
displacement  thickness.  This  displacement  effect  causes  changes  in  the 
properties  at  the  edge  of  the  boundary  layer. 
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External  vorticity  effect.  At  supersonic  and  hypersonic  speeds, 
a  blunt  body  produces  a  curved  detached  bow  shock  wave  which  in  turn 
produces  an  inviscid  rotational  flow  field.  A  boundary  layer  growing 
on  a  body  under  such  a  flow  field  would  have  different  "edge"  conditions 
because  of  the  shock-generated  inviscid  external  vorticity.  This  effect 
can  be  further  subdivided  into  an  entropy  gradient  effect  and  a  stag¬ 
nation  enthalpy  gradient  effect.  For  most  practical  problems,  it  is 
satisfactory  to  treat  the  inviscid  flow  field  as  isoenergetic,  and  thus 
under  these  conditions  only  the  entropy  gradient  effect  needs  to  be  con¬ 
sidered. 

Slip  and  temperature  jump  effects.  Under  highly  rarefied  flow 
conditions,  the  gas  tangential  velocity  at  the  wall  is  not  zero  relative 
to  the  wall  and  the  temperature  of  the  gas  at  the  wall  is  not  the  wall 
temperature  value.  A  continuum  analysis  based  on  the  Navier-Stokes 
equations  can  be  applied  if  the  wall  boundary  conditions  on  the  tan¬ 
gential  velocity  and  static  temperature  are  modified  to  include  the 
noncontinuum  effects  of  velocity  slip  and  temperature  jump. 

The  roles  of  first-  and  second-order  boundary- layer  effects 
will  be  made  clearer  in  the  following  discussion. 

II.  INTRODUCTION 

As  a  high-speed  vehicle  enters  the  earth's  atmosphere,  the  flow 
field  about  the  vehicle  changes  from  near  free-molecule  flow  in  the 
outer  regions  of  the  atmosphere  through  various  noncontinuum  regions 
to  a  more  classical  continuum  inviscid  outer  flow  and  inner  viscous 
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laminar  and  turbulent  flow  fields.  At  high  speed  the  gas  is  not  a 
simple  ideal  gas  since  chemical  equilibrium  and  nonequilibrium  effects 
on  the  thermodynamic  properties  are  also  important. 

The  part  of  this  complex  problem  which  will  be  of  concern  in 
this  dissertation  is  the  extension  of  classical  laminar  boundary- layer 
flow  theory  into  the  lower  Reynolds  number  regime.  In  modem  hyper¬ 
sonic  wind  tunnels,  the  aim  is  toward  flight  simulation  since  flight 
duplication  is  almost  never  achieved.  It  is  also  worth  noting  that  in 
these  hypersonic  test  facilities,  the  flow  fields  over  blunt-nosed 
bodies  are  almost  always  laminar  since  the  effect  of  the  blunt  nose 
reduces  the  local  Reynolds  number  in  the  flow  field,  and  thus  natural 
transition  from  laminar  to  turbulent  flow  seldom  occurs. 

The  incentive  for  this  research  began  in  1959  with  some  un¬ 
expectedly  high  drag  measurements  found  in  testing  a  slender  blunt 
cone  in  the  hotshot  tunnels  of  the  von  Karman  Gas  Dynamics  Facility 
(VKF) .  The  application  of  simple  boundary- layer  theory  failed  to  pre¬ 
dict  the  experimentally  observed  data.  Whitfield  and  Griffith  [l]1 
developed  an  approximate  theory  for  predicting  total  rero-lift  drag 
based  upon  the  boundary- layer  theory  of  Cohen  and  Reshotko  [2,3],  dis¬ 
placement  effects  based  on  the  work  of  Lees  and  Probstein  [4]  and 
and  Probstein  [5] ,  and  transverse  curvature  [TVC)  considerations  as 
treated  by  Probstein  and  Elliott  [6],  The  method  briefly  was  as 
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references  in  the  bibliography. 


5 


A  E  DC-T  R-68-148 


follows:  First,  the  friction  drag  and  displacement  thickness  were  com* 
puted  by  the  Cohen  and  Reshotko  integral  method  [3]  and  a  correlation 
of  experimental  pressure  distributions.  Then,  with  the  aid  of  the 
above-mentioned  theories  [4,5,6]  the  displacement-induced  pressure  and 
friction  drags  and  transverse-curvature-induced  friction  drag  were 
computed.  Whitfield  and  Griffith  assumed  that  the  theories  of 
Probstein  [S]  and  Probstein  and  Elliott  [6]  which  were  developed  for 
sharp  cones  could  also  be  applied  to  spherically  blunted  cones  by 
assuming  that  the  local,  induced  skin  friction  at  the  sphere-cone 
junction  was  the  same  as  would  exist  at  the  same  point  on  a  sharp  cone 
(i.e.,  tw  (blunt  sphere-cone  junction,  0C)  =  tw  (sharp,  ©c)).  In 
addition,  constant  boundary- layer  edge  properties  were  assumed  over  the 
conical  afterbody.  The  justification  for  the  noted  assumptions 
a  posteriori  was  that  the  induced  effects  were  small,  the  corrections 
were  in  the  proper  direction,  and  the  agreement  between  analysis  and 
experiment  was  good.  The  calculations  were  made  using  both  perfect 
(Y  a  1.4)  and  equilibrium  gas  models.  The  highly  approximate  nature 
of  the  theoretical  model  used  by  Whitfield  and  Griffith  [1]  was 
recognized;  however,  a  correlation  parameter  was  defined  and  a  method 
developed  which  correlated  well  the  experimental  drag  data  at  -  19. 
The  agreement  between  Whitfield  and  Griffith's  prediction  and  experi¬ 
mental  data  can  be  seen  from  Figure  1. 

The  good  agreement  which  Whitfield  and  Griffith  found  between 
their  predictions  and  the  experimental  data  was  considered  somewhat 
fortuitous  since  several  boundary- layer  effects  were  not  considered 
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Fig.  1  Comparison  of  Hypersonic,  Cold  Wall,  Zero-Lift  Drag  of  Blunt,  Slender  Cones 
in  Continuum  and  Rarefied  Flow  Regimes  (from  Whitfield  and  Griffith  [1]) 
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by  them  and  some  of  those  treated  were  considered  in  a  highly  approxi¬ 
mate  manner.  The  next  attempt  at  improving  the  theoretical  model  and 
numerical  solutions  was  made  by  Lewis  and  Whitfield  [7],  The  basic 
model  used  by  Lewis  and  Whitfield  was  the  nonsimilar  laminar  boundary- 
layer  theory  of  Clutter  and  Smith  [8]  including  approximate  transverse 
curvature  terms  and  an  inviscid  outer  flow  field.  The  viscous  boundary- 
layer  and  inviscid  outer  flow  fields  were  coupled  in  the  following 
manner.  The  inviscid  blunt  body  and  characteristics  solutions  were 
obtained  using  the  method  of  Inouye,  Rakich  and  Lomax  [9].  The  surface 
pressure  distribution  over  the  geometric  body  was  then  used  to  determine 
the  conditions  at  the  edge  of  the  boundary- layer.  The  "effective  body" 
as  seen  by  the  outer  flow  was  obtained  by  perturbing  the  geometric  body 
by  the  displacement  thickness  as  follows: 

r  --  «  r  +5*  cos  a 
eff  w 

The  inviscid  blunt  body  and  characteristics  solutions  were  obtained 
over  an  approximate  effective  body  and  another  boundary- layer  solution 
was  obtained,  This  process  was  continued  until  there  was  negligible 
change  in  the  pressure  distribution  over  the  body. 

At  the  time  of  the  calculations  by  Lewis  and  Whitfield,  they 
found  it  possible  only  to  treat  spherical-nosed  bodies  in  the  inviscid 
blunt  body  solution,  and  thus  the  spherical  nose  was  expanded  or 
contracted  by  the  displacement  thickness  at  the  starting  line  (M^  = 

1.05)  for  the  characteristics  solution.  Moreover,  approximate  curve 
fit  procedures  were  developed  for  describing  the  effective  body 
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downstream  of  the  starting  line. 

The  theoretical  model  and  numerical  treatment  used  by  Lewis  and 
Whitfield  [7]  were  superior  to  those  used  by  Whitfield  and  Griffith 
[1] ;  however,  the  results  obtained  for  prediction  of  zero-lift  drag  were 
in  poorer  agreement  with  the  experimental  data  as  shown  in  Figure  2.  A 
comparison  of  drag  components  computed  by  Whitfield  and  Griffith  and  by 
Lewis  and  Whitfield  is  shown  in  Figure  3.  The  skin-friction  drag  and 
transverse-curvature-induced  drag  components  predicted  by  Whitfield  and 
Griffith  were  each  about  15  per  cent  larger  than  those  predicted  by 
Lewis  and  Whitfield.  The  minor  drag  components  due  to  displacement- 
induced  friction  and  pressure  drags  were  in  better  agreement;  however, 
their  contributions  were  about  one  order  of  magnitude  less  than  the 
first  two  components  discussed. 

The  defects  in  both  the  analyses  by  Whitfield  and  Griffith  and 
by  Lewis  and  Whitfield  were  recognized  and  in  fact  sustained  the 
interest  for  the  work  described  herein.  Both  analyses  were  approximate 
but  based  upon  what  were  considered  some  of  the  best  theoretical  models 
and  numerical  methods  available  at  the  time  the  studies  were  made. 

Both  theoretical  models  approximately  treated  TVC  and  displacement; 

Lewis  and  Whitfield  included  the  nonsimilar  boundary- layer  effects, 
whereas  Whitfield  and  Griffith  used  the  similar  boundary-layer  solutions 
and  integral  method  of  Cohen  and  Reshotko  [2,3]. 

Those  involved  in  the  above-described  analyses  of  the  experi¬ 
mental  drag  data  were  aware  of  the  hazard  of  application  of  approximate 
analyses  outside  of  ranges  where  they  were  developed  and  tested. 
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Experimental  Data 
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AEDC-VKF  Tunnel 

T0  x  10"3,  °K 

Tw^o 

■ 

16-in.  (1) 

~17 

3. 5-4. 5 

~0. 1 

♦ 

50-in.  (H) 

16-21 

~3 

~0.1 

• 

100-in.  (F) 

~21 

~4 

-0.075 

♦ 

Low  Density  (L) 

8. 8-9. 6 

2-3 

0.20-0.25 

A 

Low  Density  (L) 

9.4 

1.67 

~0.3 

Fig.  2  Comparison  of  Experimental  and  Predicted  Zero-Lift  Drag  for  a  9-degree 
Sphere-Cone  at  M„,  =  9  and  18  (from  Lewis  and  Whitfield  [7]) 
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Fig.  3  Comparison  of  Drag  Components  According  to  Lewis  and  Whitfield's 
Analysis  with  the  Previous  Analysis  of  Whitfield  and  Griffith  (from 
Lewis  and  Whitfield  [7]) 
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The  work  of  Lewis  and  Whitfield  has  been  extended  to  include  the 
effects  of  vorticity  interaction  as  first  suggested  by  Perri  and  Libby 
[10]  and  more  fully  described  by  Hayes  [11]  and  Hayes  and  Probstein  [12]. 
The  effect  of  inviscid  external  vorticity  is  simple  conceptually.  When 
the  Reynolds  number  is  sufficiently  low,  the  laminar  boundary  layer 
entrains  the  high  temperature,  low  density  gas  which  crossed  the  curved 
bow  shock.  The  velocity  at  the  "edge"  of  the  boundary  layer  increases 
above  that  along  the  inviscid  wall  streamline,  and  the  velocity  gradient 
at  the  edge  of  the  boundary  layer  is  first  increased  and  then  decreased 
depending  on  the  length  of  the  body  and  Reynolds  number  or  boundary 
layer-to-shock  layer  thickness. 

To  help  clarify  the  physical  aspects  of  shock-generated  inviscid 
external  vorticity,  some  typical  results  from  a  condition  where 
vorticity  effects  are  large  are  presented  and  discussed.  The  uncoupled 
invisexd  and  viscous  flow  fields  over  a  15-degree  half-angle  spherically 
blunted  cone  at  ■  18  are  shown  in  Figure  4.  Lines  of  minimum  entropy 
and  "edge"  of  the  entropy  layer  are  shown  in  the  inviscid  shock  layer. 
The  term  "edge"  of  the  entropy  layer  as  used  herein  denotes  a  line  above 
which  the  entropy  is  nearly  constant  and  within  one  per  cent  of  the 
constant  conical  limit.  The  shock  wave  has  an  inflection  point  where 
the  shock  angle  is  a  minimum,  the  shock  is  weakest,  and  the  entropy 
increase  across  the  bow  shock  wave  is  a  minimum.  Since  the  entropy  is 
constant  along  streamlines,  the  entropy  will  be  a  minimum  along  the 
streamline  which  passes  through  the  shock-wave  inflection  point.  The 
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Fig.  4  Unperturbed  Inviscid  Shock  Layer  and  Viscous  Boundary  Layer  over  a 
Spherically  Blunted  15-degree  Half-Angle  Cone  at  -  18  and 
Tw/T„  =  0.066 
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classical  or  first-order  boundary- layer  thickness  is  shown  where  the 
Reynolds  number  was  chosen  such  that  the  boundary- layer  thickness 
would  be  less  than  the  shock  layer  thickness,  and  the  boundary- layer 
displacement  thickness  is  also  shown.  Although  not  distinguishable  in 
this  figure,  the  boundary- layer  displacement  thickness  is  negative 
over  the  spherical  nose.  For  the  conditions  shown  in  Figure  4,  the 
boundary  layer  is  thinner  than  the  shock  layer  over  the  entire  body 
and  thinner  than  the  entropy  layer  in  the  nose  region,  but  the  boundary 
layer  becomes  thicker  than  the  entropy  layer  for  s  >  26.  Downstream 
of  this  point  the  thickness  of  the  entropy  layer  decreases  monotonically 
whereas  the  shock  and  boundary-layer  thicknesses  increase.  It  should 
be  emphasized  that  the  flow  fields  shown  in  Figure  4  were  calculated 
on  the  basis  of  being  completely  uncoupled,  and  the  conditions  were 
chosen  to  illustrate  some  of  the  features  of  the  flow  fields. 

The  results  from  the  characteristics  solution  were  interpolated 
normal  to  the  body  at  several  axial  stations,  and  the  resulting 
velocity,  temperature,  pressure,  and  entropy  distributions  are  shown 
in  Figures  5-8.  The  results  were  nondimensionalized  by  the  value  of 
the  variables  along  the  inviscid  wall  for  convenience  of  presentation. 
Since  the  sphere-cone  tangent  point  was  at  z  =  0.741,  the  profiles  at 
z  =  0.75  represent  supersonic  flow  (M^  ■  2,15  at  z  ■  0.75)  over  the 
spherical  nose  where  the  effects  of  longitudinal  curvature  on  the 
normal  pressure  gradient  are  evident.  The  flow  field  near  the  wall 
adjusts  quickly,  but  the  nose  influence  persists  until  z  >  7.5  in  the 
flow  field  near  the  shock  wave.  For  thin  boundary  layers  over  the 
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Fig.  8  Temperature  Distribution  Through  an  Inviscid  Unperturbed  Shock  Layer  Over 
a  15-degree  Half-Angle  Cone  at  M*,  =  18 
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conical  afterbody  (fi/ySh  <  0.3),  the  usual  boundary- layer  assumption  of 
3p/3y  =  0  is  supported  by  these  data. 

Primary  interest,  however,  in  these  data  is  in  the  development 
of  the  velocity  and  temperature  profiles  and  their  relationship  to  the 
entropy  profiles.  The  formation  of  the  entropy  layer  can  be  seen  from 
Figure  6.  The  monotonic  increase  in  the  velocity  and  decrease  in  the 
temperature  gradients  at  the  wall  can  be  seen  from  Figures  7  and  8. 

The  effects  of  these  gradients  in  second-order  boundary-layer  theory 
will  be  discussed  in  Chapter  II.  Far  downstream  on  the  cone  (z  >  20) 
and  even  for  thin  boundary  layers  (6/ysh  <  0.5),  the  inviscid  velocity 
at  the  "edge"  of  the  boundary  layer  would  be  about  30  per  cent  greater 
than  the  inviscid  wall  value,  and  the  temperature  would  only  be  about 
20  per  cent  of  the  inviscid  wall  value.  These  large  changes  in  edge 
conditions  are  not  usually  considered  in  boundary- layer  analyses. 

It  is  important  to  note  that  the  15-degree  half-angle  cone  was 
chosen  to  magnify  the  effects  of  the  entropy  layer  by  compressing  the 
effects  nearer  the  nose  than  would  be  the  case  for  smaller  half-angle 
cones.  However,  in  addition  to  considering  these  effects  on  the 
inviscid  flow  field,  it  is  important  to  consider  the  effects  of  the 
coupling  between  the  entropy  and  boundary  layers. 

The  coupling  of  the  inviscid  outer  flow  and  viscous  boundary 
layer  flow  complicates  the  picture  described  above,  and  some  of  the 
difficulties  caused  by  this  coupling  are  considered  in  this  disser¬ 
tation.  Also,  at  sufficiently  low  Reynolds  numbers  noncontinuum  effects 
should  be  considered  especially  as  they  affect  the  wall  boundary 
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conditions  on  the  velo'city  and  temperature  profiles. 

About  the  time  of  the  presentation  of  the  results  by  Lewis  and 
Whitfield  [7],  Davis  and  Flugge-Lotz  [13]  published  a  report  on  second- 
order  boundary- layer  effects  on  hyperboloids  and  paraboloids  at 
infinite  Mach  number  and  spheres  at  a  Mach  number  of  ten.  The  theory 
of  Van  Dyke  [14]  was  used  with  an  implicit  finite  difference  scheme 
originally  proposed  by  Flugge-Lotz  and  Blottner  [15]  for  treating  the 
classical  first-order  boundary  layer  over  two-dimensional  bodies. 

As  will  be  shown  in  this  dissertation,  the  theory  of  Van  Dyke  when 
coupled  with  the  implicit  finite  difference  method  of  Davis  and 
Flugge-Lotz  gives  a  powerful  tool  for  extending  classical  boundary- 
layer  theory  to  lower  Reynolds  number. 

More  recently  Adams  [16]  used  the  Levy-Lees  transformation  [17] 
to  solve  the  first-  and  second-order  boundary- layer  theory  again  with 
the  finite  difference  method  of  Davis  and  Flugge-Lotz,  Adams  also 
developed  an  approximate  method  based  on  the  Newtonian  theory  for  the 
pressure  distribution  over  hyperboloids  and  paraboloids.  The  essential 
difference  between  the  second-order  treatments  of  Davis  and  Flugge-Lotz 
and  that  of  Adams  was  in  the  treatment  of  the  displacement  effect. 

Adams  was  able  to  treat  analytic  bodies  of  arbitrary  length  whereas 
Davis  and  Flugge-Lotz  were  restricted  to  the  nose  region  since  in  their 
method  the  effective  body  was  approximated  by  shifting  and  expanding 
the  original  body  by  the  displacement  thickness.  Although  the  shifted 
and  expanded  body  technique  was  useful  and  had  some  advantages  over 
other  methods,  it  was  limited  to  the  nose  region  since  it  did  not  provide 
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satisfactory  fits  of  the  effective  body  over  regions  far  removed  from 
the  stagnation  region.  The  method  developed  by  Adams  was  limited  to 
analytic  bodies  since  it  was  based  on  the  Newtonian  impact  theory  and 
was  thus  not  applicable  to  spherically  blunted  cones. 

A  method  has  recently  been  developed  by  Marchand,  Lewis,  and 
Davis  [18]  for  treating  the  second-order  displacement  effect  over 
nonanalytic  shapes  such  as  spherically  blunted  cones.  The  method  used 
by  Marchand,  Lewis,  and  Davis  for  treating  the  displacement  effect 
will  be  described  later  herein,  and  the  remaining  second-order  effects 
followed  Davis  and  Flugge-Lotz  exactly  except  for  the  effects  of 
vorticity  which  were  coupled  to  the  displacement  effect. 

The  primary  goal  of  this  dissertation  will  be  to  present  a  com¬ 
parison  of  a  first-order  treatment  of  second-order  boundary- layer 
effects  with  treatments  of  Van  Dyke's  second-order  theory  and  a  com¬ 
parison  of  both  treatments  with  experimental  data.  From  these  com¬ 
parisons  some  conclusions  will  be  drawn  regarding  the  ranges  of  appli¬ 
cability  of  the  theories. 

In  Chapter  II  the  theory  of  higher-order  boundary- layer  effects 
is  given.  A  review  of  Van  Dyke's  second-order  theory  is  presented 
including  a  discussion  of  the  asymptotic  matching  principle.  The  source 
of  the  various  second-order  effects  is  indicated,  and  the  coupling 
between  second-order  vorticity  and  displacement  is  shown. 

The  first-order  treatment  of  higher-order  boundary- layer  effects 
is  given  in  Chapter  III.  The  effects  of  vorticity,  slip  and  tempera¬ 
ture  jump  [considered  as  a  single  effect),  and  displacement  are  described. 
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The  effect  of  first-order  transverse  curvature  was  developed  by  Clutter 
and  Smith,  and  the  approximate  nature  of  the  treatment  caused  by  an 
error  in  the  original  theoretical  development  is  discussed. 

Numerical  results  and  discussion  will  be  presented  in  Chapter 
IV.  The  conditions  treated  by  the  first-  and  second-order  theories 
will  be  given,  and  the  numerical  results  will  be  presented  and  dis¬ 
cussed. 

Comparisons  of  the  numerical  results  from  the  first-  and 
second-order  theories  will  be  made  with  experimental  data  and  the 
results  discussed  in  Chapter  V. 

The  conclusions  from  the  study  will  be  given  in  Chapter  VI. 

From  the  comparison  of  numerical  results  and  experimental  data,  the 
ranges  of  applicability  of  the  theories  will  be  discussed,  and  the  need 
for  future  work  will  be  indicated. 

Separate  effects  and  problems  will  be  treated  in  appendixes.  The 
numerical  treatment  of  the  first-order  boundary- layer  theory  is  des¬ 
cribed  in  Appendix  A.  The  method  used  to  treat  first-  and  second-order 
displacement  effects  is  described  in  Appendix  B.  Finally,  the  effect 
of  vorticity  on  the  displacement  thickness  is  derived  in  Appendix  C. 
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CHAPTER  II 

REVIEWOF  SECOND-ORDER  BOUNDARY-LAYER  THEORY 


The  second-order  boundary- layer  theory  of  Van  Dyke  [14]  is 
reviewed  to  indicate  the  salient  features  of  the  theory  and  the  roles 
of  the  various  second-order  effects.  The  critical  role  of  the  match¬ 
ing  of  the  inner  and  outer  expansions,, based  on  Van  Dyke’s  asymptotic 
matching  principle,  is  discussed. 

In  the  second-order  boundary-layer  theory,  Van  Dyke  [14] 
identified  seven  separate  second-order  effects  which  were  longitudinal 
curvature,  transverse  curvature,  slip,  temperature  jump,  entropy 
gradient,  stagnation  enthalpy  gradient,  and  displacement.*  The 
separation  of  the  effects  of  vorticity  and  displacement  is  somewhat 
arbitrary;  however,  the  presentation  given  below  follows  Van  Dyke's 
treatment . 

I.  SECOND-ORDER  BOUNDARY-LAYER  EQUATIONS 

The  Navier-Stokes  equations  are  written  in  the  orthogonal 
coordinate  system  (s,n)  shown  in  Figure  9,  where  s  is  the  non- 
dimensional  distance  along  the  surface,  n  is  the  nondimensional  distance 
normal  to  the  surface,  <(s)  is  the  longitudinal  surface  curvature,  and 


brief  physical  explanation  of  each  of  the  second-order 
boundary- layer  effects  is  given  in  Chapter  I. 
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rw(s)  is  the  radius  of  the  axisymmetric  body  where  all  lengths  are  non- 
dimensionalized  by  the  nose  radius  of  curvature  a*.  Further,  all  veloci¬ 
ties  are  referred  to  U*,  pressure  to  P*U*2,  density  to  p*,  temperature 
to  U*2/c*,  entropy  to  c*,  enthalpy  to  U*2,  and  viscosity  to  the  value 
of  v  at  T*  =  U*2/c*. 

00  p 

The  continuity  equation  becomes 

(rjpu)s  +  (krJpv)n  =0  (1) 

where  subscripts  denote  differentiation  and  j  =>  0  for  two-dimensional 
flow  and  j  ■  1  for  axisymmetric  flow.  The  radius  r  =  rw  +  n  cos  a  and 
the  function  k  =  1  +  <n  are  used  for  convenience. 

The  Navier-Stokes  momentum  equations  become 

e'2  [p(uus/k  +  vi^  +  Kuv/k)  +  ps/k]  ■  {y[un  +  (vs  -  Ku)/k])n 

+  (2/k)  [u(us  +  KV]/k]s  +  y(2«/k  +  j  cos  a/r) [u„  +  (vs  -  *u)/k] 

+  [2jyrs/(kr)]  [ (us  *  <v)/k  -  urs/Ckr)  -  v  cos  a/r] 

+  (1/k]  U[(us  +  i«v)/k  +  vn  +  j  Curs/k  +  v  cos  a)/r]}s  (2) 

and 

e-2  [p (uvs/k  +  wn  -  ieu2/k)  +  pn]  *  2(yvn)n  +  {y[Un  +  (vs  -  Ku)/k])s 
+  2y(ic/k  +  j  cos  a/r)vn  -  2uk(us  +  icv)/k2 
-  2jy (urs/k  +  v  cos  a]/r2  +  jyrs  [un  +  (vs  -  Ku)/k]/{kr) 

+  U[(us  +  <v)/k  +  vn  +  j  (urs/k  +  v  cos  a)/r])n  (3) 
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The  energy  equation  is 

e"2  [P(uTs/k  +  vTn)  -  (ups/k  +  vpn)]  =  o"1  (pTs/k)s/k 
+  o"1  (uTn)n  +  jo"1  uTsrs/(k2r)  +  o-1  (ic/k  +  j  cos  a/r)  Tn  +  *  (4) 

where  a  is  the  Prandtl  number,  and  the  dissipation  function  is 
*  =  jj{2[(us  +  Kv)/k]2  +  2vn2  +  2jr"2  Curs/k  +  v  cos  a)2 
+  [Un  +  (vs  -  Ku)/k]2} 

+  *  [Cus  +  *v)/k  +  vn  +  jr-1  (urg/k  +  v  cos  a)]2  (5) 

The  nondimens ional  equation  of  state  is 

p  =  (Y-l)pT/y  (6) 

The  viscosity  law  is 

y  =  y(T)  (7) 

The  second  coefficient  of  viscosity  is  A. 

Van  Dyke  defined  the  expansion  parameter  £  to  be 


If  a  power  viscosity  law  y  «  Tw  is  assumed,  then 
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Thus  with  y  »  1.4,  e  =  0.632S  M^/Re.,^/2  for  the  linear  viscosity  law 
(u  *  1.0),  and  e  =  0.795  (M^/Re^) V2  for  the  square  root  viscosity  law 
(to  *  0.5). 2 

The  outer  (outside  the  boundary  layer)  flow  variables  in  the 
Navier-Stokes  equations  (including  continuity  and  energy  equations  in 
the  system)  are  expanded  in  powers  of  e  as  follows: 


u(s,n;e)  U1(s,n)  ♦  eU2(s,n)  ♦  .  .  .  (9a) 
v(s,n;e)  v  V^s.n)  ♦  eV2(s,n)  ♦  .  .  .  (9b) 
p(s,n;e)  •v  P1(s,n)  ♦  eP2(s,n)  ♦  .  .  .  (9c) 
p(s,n;e)  v  Rj(s#n)  +  eR2(s,n)  +  .  .  .  (9d) 
T(s,n;c)  v  T^Sjn)  +  eT2(s,n)  ♦  .  .  .  (9e) 


Substituting  Equations  9  into  Equations  1-6  and  equating  terms 
in  e°  gives 

(rJRjU^g  ♦  (kr^Vj),,  -  0  (10) 
Rl(UiUis/k  ♦  VjUm  ♦  KUjVj/k)  ♦  Pl5/k  =»  0  (11) 
Rl(UiVls/k  ♦  VjVln  -  xU^/k)  ♦  pln  -  0  (12) 


2With  linear  viscosity  law,  the  expansion  parameter_is  related 
to  the  hypersonic  interaction  parameter  vw,  by  *  *  0.6325  v,,,  x 


27 


AEDC-TR.68-148 


RlOJjTig/k  ♦  VjT^}  -  (UjPls/k  ♦  VjPln)  =  0  (13) 

and 


Pj  =  (y-1)  RjTj/y 


(14) 


which  are  the  first-order  inviscid  flow  equations. 

The  stream  function  is  defined  in  the  usual  way 


and 


*JR1U1  *  *ln 


krJRlVl  *  ‘  *ls 


(15) 


The  entropy  and  stagnation  enthalpy  are  constant  along 
streamlines  and  can  be  written 

2 

sx  =  SjC rT)f  .  and  Hj  -  Tx  ♦  IuA  =  H^)  (16) 

The  second-order  equations  are  obtained  by  equating  terms  in 
e1,  and  the  tangential  momentum  equation  is 


R1lCulU2)s/k 


VlU2n 


*  V2Uln 


<  (uxv2  ♦ 


u2v1)/k] 


+  R2(UlUls/k  +  VlUln  +  KUjVj/k)  +  P2s/k  =  0 


(17) 


This  equation,  typical  of  the  system  of  equations,  is  linear  in 

the  second-order  terms  and  inviscid  since  viscous  terms  in  the  Navier- 

2 

Stokes  equations  are  of  third-order  being  multiplied  by  e  .  Thus,  a 
third-order  expansion  would  introduce  viscous  terms  in  the  equations 
for  the  outer  flow. 
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The  normal  derivatives  of  Vj  and  Pj  at  the  surface  are  obtained 
by  evaluating  Equations  10  and  12  at  n  =  0,  and  one  finds 

vl„  *  -  tCrA)'1  (r„jHlU,)sJn=0  (18) 

and 

Pin  *  '  <19> 

The  normal  component  of  the  tangential  velocity  Uj  is  obtained 
from  the  vorticity  of  the  outer  flow  which  is 

fil  =  '  Uln  +  (Vis  *  KUl)/k  (2°) 

Crocco's  vortex  theroem  [12]  is 

x  iij  =  Vllj  -  TjVSj  (21) 

at  the  surface,  and  substituting  Equations  15  and  16  into  Equation  21 
gives 

nl  *  rjRi  C^Sj'-  Hj»)  at  n  -  0  (22) 

where  prime  denotes  d/dVj. 

Substituting  Equation  20  into  Equation  22,  one  can  obtain 

uln  =  -  kUj  -  v’RjP’jSj'-  Hj')  (22a) 

to  be  evaluated  at  n  =  0. 

The  normal  gradients  of  Tj  and  Ri  are  obtained  by  differentiat¬ 
ing  Equations  16  and  14  with  respect  to  n,  using  Equation  22a,  and 
evaluating  the  results  at  the  wall  yields 
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Tln  =  [<ul  *  V  WiV^O 


(23) 


and 


K  U1  -j  ,\ 

Rln  * 

Rlul 

7^17  -  V  Risi' 

(24) 


The  outer  expansion  fails  at  the  surface  since  the  highest 
derivatives  in  the  Navier-Stokes  equations  were  lost  in  the  expansion. 
Therefore,  an  inner  expansion  is  needed  which  is  valid  near  the  surface 
although  it  will  be  invalid  far  from  the  wall.  The  inner  expansion 
is  obtained  by  defining  an  expansion  of  the  form 


u(s,n;c) 

-v  Ujts.n)  + 

eu2(s,N)  ♦  .  .  . 

(25) 

v(s.n;e) 

"V  eV^S.N)  + 

e2v2(s,N)  +  .  .  . 

(26) 

p(s,n;e) 

^  Pj (s,N)  «■ 

cp2(s,N)  ♦  .  .  . 

(27) 

P (s,n;e) 

*  Pj(s,N)  + 

ep2 (s,N)  ♦  .  .  . 

(28) 

T(s,n; e) 

*  tj (s,N)  + 

ct2(s,N)  +  .  .  . 

(29) 

where  n  =  eN. 

The  viscosity  is  expanded  in  a  series  in  T  as  follows: 

u(T)  *v  y(tj)  +  eu'(tj)t2  ♦  .  .  .  (30) 

Substituting  the  inner  expansion  Equations  25-30  into  the  Navier- 
Stokes  Equations  1-7  and  equating  terms  in  e°  gives  the  first-order 
equations,  which  are 
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Crwjpiuj)s  +  (rJpjV^N  -  0  (31) 

Pl(uiuls  4  Vlu1N)  ♦  p1$  -  Cpu1n)n  =  0  (32) 

P1N  a  0  (33) 

MVlS  +  VlN^  “  (ulPls  *  VlPlN^ 

“  o  *  ^U1N  a  ®  (34) 

and 

P,-  P^,  (35) 

which  are  Prandtl's  compressible  boundary- layer  equations. 

Equating  terms  in  e1  yields  the  second-order  boundary- layer 
equations,  which  are 

tV>  (Plu2  +  P2U!)]S  ♦  [rwj  (P jv 2  +  P2V1>]N 

*  ^ort  +  CLC  +  CTVC  *  Cdisp  +  CSTJ  <36) 

pl<ulu2s  +  u2uls  +  V1U2N  +  V2U1N}  *  p2Cululs  +  V1U1N> 

-  (uu2N  ♦  u'u1Nt2)N 

s  NWt  +  mLC  +  mTVC  +  Mdisp  +  MSTJ  t37^ 

pl(ult2s  +  u2tls  *  vlt2N  +  v2tlN^  +  ^2^1*18  +  vltlN^ 

-  u2pis(s,°)  -  (p°‘lt2)NN  -  2wu1Nu2N  "  u’ulN2t2 

*  Evort  +  ELC  +  eTVC  +  Edisp  +  ESTJ  C38) 
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P2  *  CPit2  ♦  P2ti3 

and 

P2N  =  <P!«!2 

where 

cvort  3  cdisp  =  CSTJ  “  0 
CLC  *  ‘  <T»  P1V1>N 

ctvc  *  "  cos“  Npiui^s  +  Crwj  l  cosot  NpiviJn1 

"vort  *  -  vWTlS1'V2]„=0 

mLC  ■  k[N^u1N>N  +  UU1N  *  P'VlN 

-  Np1vJu1N  -  PjUjVjJ 

-  AlJs 

^TVC  =  i  cos  (I  V  u1N/rw 
Mdisp  =  r*l(ulu2)S  +  Wls^O 
mstj  a  0 

Evort  =  rwjullRl  Tlsl'v2ln=0 

eLC  =  K[NpiUitls  ♦  o~lv  t1N  ♦  Plu12v1 

-  NuiPis(s,0)  -  2y  Uju1N] 

+  Uj [k ^i]s 


(39) 

(40) 

(38a) 

(38b) 

(38c) 

(39a) 


(39b) 

(39c) 

(39d) 

(39e) 

(40a) 


(40b) 
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H-rve  *  cos  a  W  T]Nj/rw  (40c) 

Edisp  "  *  “ltRlCUiU2)s  *  RaUiULJn.o  (40d) 

Estj  *  0  (40e) 

and 

DO 

Aa  =  N  +  /  [1  -  Pi^j2/ CR,Uj2)n_0]dN.  C 4 Of ) 

N 

The  subscripts  vort,  disp,  LC,  TVC,  and  STJ  denote  contributions 
due  to  vorticity,  displacement,  longitudinal  curvature,  transverse 
curvature,  and  slip  and  temperature  jump,  respectively.  The  differ¬ 
ential  equations  are  linear,  but  not  all  of  the  effects  are  independent 
since,  for  example  as  will  be  discussed  below,  the  V2(s,0)  in  liquation 
30a  is  related  to  the  displacement  effect  through  the  displacement 
thickness  6*. 

The  second  coefficient  of  viscosity  A  is  of  third-order  (that 
is,  in  terms  multiplied  by  f  )  and  thus  does  not  appear. 

Neglecting  surface  mass  transfer,  the  boundary  condition  along 
the  surface  requires 


v(s,0)  =  0 


(41) 


Van  Dyke  [14]  obtained  the  remaining  surface  boundary  conditions  from 
Street  [20]  including  first-order  slip-flow  conditions  as  follows: 


u(s,0)  = 


xlL 

Y 


T 


1/2 


‘‘n 


+  b, 


Y~l. 

Y 


(42) 


n=0 


and 
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T(s,0)  »  Tw  + 


(43) 


where  *  (ir/2)*/^,  bj  =  3/4,  and  c^  =  (15/8)  (ir/2)^/2  when  a  hard 
sphere,  monatomic  gas  model  with  unit  accommodation  coefficients  is 
assumed . 

If  the  surface  temperature  is  prescribed  as  Tw(s),  then  sub¬ 
stitution  of  the  inner  expansion  Equations  24-30  into  the  wall  boundary 
conditions  Equations  41-43  gives 


uj (s ,0)  *  Vjfs.O)  =  0  (44) 

u2(s,0)  -  |“  tjJ  1  u1n]n=o  (45) 

v2(s,0)  *  0  (46) 

tj(s.O)  =  tw(s)  ‘  (47) 

t2C=.°)  *c,[^  ti]1/,tuj(|>#  C48, 


Van  Dyke  [14]  gives  corresponding  expressions  when  the  tempera¬ 
ture  gradient  at  the  wall  is  prescribed. 

II.  ASYMPTOTIC  MATCHING  CONDITIONS 

The  outer  boundary  conditions  will  depend  on  the  problem  to  be 
solved;  however,  the  outer  expansion  fails  to  satisfy  the  wall  boundary 
condition  and,  in  general,  the  inner  expansion  fails  as  N  ->».  It  is 
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fundamental  to  the  method  of  matched  asymptotic  expansions  that  there 
exist  an  overlap  region  where  both  inner  and  outer  expansions  are 
valid,  and  thus  both  expansions  can  be  matched  in  this  region. 

Van  Dyke  [19]  described  the  asymptotic  matching  principle  as 
follows: 

The  m-term  inner  expansion  of  (the  n-term  outer 
expansion)  =  the  n-term  outer  expansion  of  (the 
m-term  inner  expansion).  (49) 

By  his  definition,  the  m-teim  inner  expansion  of  the  n-term 
outer  expansion  is  obtained  by  rewriting  n-terms  of  the  outer  expan¬ 
sion  in  inner  variables  (N  in  second-order  boundary- layer  theory), 
expanding  asymptotically  in  small  e,  and  truncating  the  result  after 
in  nonzero  terms.  The  converse  of  the  procedure  is  used  on  the  right- 
hand  side  of  equation  49. 

To  demonstrate  the  matching  procedure,  the  tangential  velocity 
component  is  considered  first.  Applying  the  matching  principle  with 
n  *  m  =  1  gives  the  matching  conditions  for  the  first-order  boundary- 
layer  theory.  According  to  the  matching  principle, 


1-term  outer  expansion: 
rewritten  in  inner  variables: 
expanded  for  small  c: 

1-term  inner  expansion: 

1-term  inner  expansion: 
rewritten  in  outer  variables: 


Uj (s  ,n) 

UjCs.eN) 

Uj(s.O)  *  NeUln(s,0)  + 

Ul(s,0) 

u1(s,N) 

A 

u1(s,n/c) 
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expanded  for  small  e:  Uj(s,'°) 

1-term  outer  expansion:  Uj(s,~) 


Equating  the  results  one  obtains 

Uj(s,N)  U^(s,0)  as  N  +  «  (50a) 

In  the  same  manner,  one  finds  the  following  matching  conditions: 

Pl(s,N)  %  P^s.O)  (50b) 

PjCs.N)  ^  RjCs.O)  >  as  N  +  ■  (50c) 

t^s.n)  a,  Tj (s,0)  (50d) 


The  behavior  of  as  N  +  «  is  found  by  solving  the  first-order 
boundary- layer  equations. 

The  matching  with  m  =  1  and  n  *  2  gives  the  matching  condition 
due  to  the  displacement  thickness.  Again  applying  the  asymptotic 
matching  principle  to  the  normal  velocity  component  v  yields  the  follow¬ 
ing: 

1- term  inner  expansion:  eVj(s,N) 

rewritten  in  outer  variables:  EV^(s,n/e) 
expanded  for  small  e:  ev^(s,«0 

2-tcrm  outer  expansion:  evj(s,<») 

2- term  outer  expansion:  V1(s,n)  ♦  eV2(s,n) 

rewritten  in  inner  variables:  Vjts.eN)  +  eV2(s,eN) 
expanded  for  small  e:  Vj(s,0)  +  eNVjn(s,0)  +  .  .  , 

+  e[V2(s,0)  +  cNV2n(s,0)  +  ...J 
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1-term  inner  expansion:  e^Cs.O)  +  NVjn(s,0)] 

Equating  the  results  gives 

V2(s,0)  =  lim  [Vl£s,N)  -  NVlnCs,0)]  (51) 

N  «® 

where  use  was  made  of  the  boundary  condition  Vj(s,0)  =  0. 

Matching  with  m  «  n  «  2  gives  the  remaining  matching  conditions 


as  follows: 

u2(s,N)  U2(s,0)  ♦  NUln(s,0)  (52a) 

p2(s,N)  -v.  P2(s,0)  ♦  NPln(s,0)  (52b) 

P2(s,N)  ^  R2(s,0)  ♦  NRln(s,0)  (52c) 

t2(s,N)  a,  T2(s,0)  +  NTln(s,0)  (52U) 

Vj  (s,N)  -v  V2(s,0)  ♦  NVln(s,0)  (52e) 


where  each  of  Equations  52  is  evaluated  as  N  •> 

Equation  51  can  be  related  to  the  boundary -layer  displacement 
thickness  6*  as  follows:  The  definition  of  the  displacement  thickness 
is 

00 

6*  =  /  [1  -  Piu1/(R1U1)  0JdN  (53) 

o 

Using  the  inner  expansion  stream  function  'I'jnCSjN)  ■  r^pjUj,  and 
tmfr.O)  =  *in(s,N)  as  N  +  »,  one  obtains 
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6*  =  lim  [N  -  ^(s.NJ/^Cs.O)]  (54) 

N  ■+•  “ 

Multiplying  by  rw^Rj(s,0)U^(s,0)  =  ^^(SiO)  and  differentiating  with 
respect  to  s  yields 

(r„3R1(8#0)U1(si0)«*)s 

-  lim  [N*laN(s,N)  -  *ls(s,N)] 

N  •*  « 

■  lim  [N(-rw3p1v1)N  +  v’pjvJ 

N  -*■  « 

=  lim  [r^pjVj  -  twjN(Piv1n  +  vlPlN)J  (55) 

N  -*■  • 

where  the  last  term  in  the  last  expression  is  zero  since  the  gradient 
of  the  first-order  density  vanishes  as  N  -►  The  final  result  is 

(rwjR1(s,0)U1(s,0)6*)s  «  rwjR1(s*0)V2(s,0)  (56) 

where  Equation  51  has  been  used.  This  final  expression  provides  a 
means  of  evaluating  the  second-order  normal  velocity  component  V2(s,0) 
once  the  first-order  outer  flow  and  displacement  thickness  are  known. 

As  described  by  Van  Dyke  [ 14 J  and  Davis  and  Flugge-Lotz  [13],  this 
relation  is  important  in  second-order  boundary- layer  theory  since  it 
relates  the  displacement  effect  through  6*  with  the  vorticity  effect 
through  V2(s,0), 

It  is  instructive  to  consider  the  matching  conditions  in  a 
slightly  different  way  which  indicates  the  second-order  effects 
involved.  From  Davis  and  Flugge-Lotz,  these  conditions  are  as 
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f ol lows : 

u2(s,N)  x  -  r|>[R1T1si]n.0  -  rillUjln.0  +  [U2]n=0 

and 

t2(s,N)  -  rjNtRjTjUjSllH.o  .  KN[U1^lnBQ  .  [T2)n.„ 

to  be  evaluated  as  N  +  »,  and  where  the  terms  on  the  right-hand  side  of 
each  expression  are  due  to  vorticity,  longitudinal  curvature  and 
displacement,  respectively. 

This  rather  extensive  review  of  second-order  boundary- layer 
theory  has  been  given  to  indicate  the  various  second-order  effects  and 
their  relationship  to  first-order  theory.  The  critical  role  of  the 
matching  of  inner  and  outer  expansions  has  also  been  indicated  in 
addition  to  the  coupling  which  exists  between  the  displacement  and 
vorticity  effects. 

The  first-  and  second-order  equations  were  solved  numerically 
by  Van  Dyke  [14]  at  the  stagnation  point  of  a  sphere  at  infinite  Mach 
number.  Later  Davis  and  Fliigge-Lotz  [13]  using  an  implicit  finite 
difference  method,  solved  the  stagnation  region  for  hyperboloids  and 
paraboloids  at  infinite  Mach  number  and  for  a  sphere  at  M,,,  ■  10.  The 
conditions  considered  by  Van  Dyke  gave  a  displacement  thickness  near 
zero,  and  he  neglected  that  effect.  The  results  of  the  method  used 
by  Davis  and  Fliigge-Lotz  for  treating  the  displacement  effect  (the 
shifted  and  expanded  body)  was  independent  of  the  expansion  parameter 
e,  and  thus  the  results  could  be  scaled  in  e.  This  result  is  generally 
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not  possible  to  obtain  because  the  effects  of  displacement  are  depend¬ 
ent  on  the  perturbation  of  the  inviscid  outer  flow  which  is  a  global 
rather  than  local  effect. 

Because  of  the  coupling  of  displacement  and  vorticity,  the  effect 
has  been  treated  as  either  a  displacement  or  vorticity  effect.  The 
effect  is  more  clearly  seen  by  considering  the  second-order  inviscid 
pressure  gradient  which  is 

P2sCs.O)  -  -[R1C(J1U2)3  4-  R2UlUls  -  r1)VT1V2s;in.O  <s9) 

where  the  last  term  arises  due  to  the  coupling  between  vorticity  and 
displacement  (see  Equation  56).  According  to  the  second-order  theory, 
the  term  may  be  considered  as  either  a  displacement  or  vorticity 
effect.  Van  Dyke  [21]  calls  the  effect  a  "displacement  pressure"  effect 
if  the  term  is  considered  as  a  displacement  effect  and  as  a  "displace¬ 
ment  speed"  effect  if  treated  as  a  vorticity  effect.  Both  Van  Dyke  [14] 
and  Davis  and  Flugge-Lotz  treated  the  term  as  a  displacement  speed 
(vorticity)  effect.  Recently  Adams  [16]  and  Marchand,  Lewis,  and  Davis 
[18]  treated  analytic  bodies  and  spherically  blunted  cones,  respectively, 
by  both  methods.  It  was  observed  first  by  inspection  of  the  numerical 
results  and  later  proved  by  Adams  [16]  that  the  displacement  pressure 
treatment  results  were  not  unique  but  were  affected  by  where  the  outer 

7 

boundary  conditions  were  imposed.  The  displacement  speed  treatment, 

3The  term  unique  will  be  used  to  denote  second-order  results 
which  were  not  affected  by  the  location  of  the  outer  edge  of  the  boundary 
layer  or  the  matching  point  for  the  inner  and  outer  flow  field  solutions. 
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however,  gave  results  independent  of  the  matching  point  or  edge  of  the 
first-order  boundary  layer.  The  disadvantage  of  the  displacement  speed 
treatment  was  that  it  gave  unrealistic,  separate  vorticity  and  dis¬ 
placement  effects.  Adams  [16]  and  Marchand,  Lewis,  and  Davis  [18] 
suggested  that  the  only  proper  physical  treatment  was  a  combined 
vorticity-displacement  interaction.  A  significant  disadvantage  in  the 
combined  vorticity-displacement  interaction  is  that  it  makes  comparison 
with  other  theoretical  or  numerical  treatments  more  difficult.  It 
should  be  noted  that  although  the  displacement  pressure  results  were 
found  to  be  nonunique,  comparison  of  results  for  the  sum  of  displace¬ 
ment  plus  vorticity  effects  between  the  two  treatments  showed  errors 
of  about  S  per  cent  in  the  displacement  pressure  results  for  the  wall 
shear  stress.  Because  the  errors  found  by  Adams  [16]  and  Marchand, 
Lewis,  and  Davis  [18]  were  small  and  since  the  displacement  speed 
treatment  gave  unrealistically  large  displacement  effects  (treated  as 
a  separate  effect] ,  the  second-order  results  used  for  comparison  with 
the  first-order  treatment  of  vorticity  and  displacement  came  from  a 
displacement  pressure  treatment. 
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CHAPTER  III 

FIRST-ORDER  TREATMENT 

The  first-order  treatment  of  higher-order  boundary- layer  effects 
presented  herein  is  based  on  a  modification  to  the  nonsimilar  boundary- 
layer  theory  of  Clutter  and  Smith  [8],  The  modification  includes 
addition  of  the  separate  or  combined  effects  of  the  vorticity  inter¬ 
action  theory  of  Hayes  and  Prob stein  [12]  and  the  slip  and  temperature 
jump  wall  boundary  conditions  of  Van  Dyke  [14].  In  addition,  the 
approximate  treatment  of  transverse  curvature  by  Clutter  and  Smith  and 
the  effects  of  displacement  are  also  considered. 

An  error  exists  in  the  transverse  curvature  term  in  the  Clutter 
and  Smith  theory  which  is  retained  in  the  present  treatment.  The 
reasons  for  not  correcting  the  error  were  twofold.  First,  to  correct 
the  error  without  redefining  the  transformation  variables  would  have 
introduced  several  additional  terms  and  substantially  complicated  the 
numerical  solutions.  The  error  has  recently  been  corrected  by  Jaffe, 
Lind,  and  Smith  [22]  by  means  of  new  transformation  variables.  The 
second  reason  for  not  changing  the  Clutter  and  Smith  variables  was  the 
widespread  use  of  the  computing  machine  program  throughout  the  country 
in  engineering  applications.  It  was  felt  desirable  to  establish  the 
error  caused  in  wall  shear  stress,  for  example,  using  the  Clutter  and 
Smith  method.  Comparisons  will  be  given  in  a  later  chapter  which  will 
show  the  magnitude  of  the  error  in  the  integrated  friction  drag. 
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The  first-order  boundary- layer  equations  are  given  in  this 
chapter  with  derivations  of  the  effects  of  vorticity  and  slip  and 
temperature  jump  (STJ).  The  method  of  solution  of  the  first-order 
boundary-layer  equations  is  described  in  Appendix  A. 

I.  BOUNDARY-LAYER  EQUATIONS 

The  first-order  steady-flow  equations  for  an  axisymmetric 
laminar  boundary  layer  in  the  coordinate  system  shown  in  Figure  10 
including  transverse  curvature  are  as  follows: 


(rpu)x  +  (rpv)y  =  0 

(60) 

p(uux  +  vuy)  =  -  px  +  H  ryuy  +  (uuy)y 

(61) 

p(uHx  -  vHy)  =  i  ry  fjL-  Hy  +  y(l  -  1/Pr)uuy 

+  [fe  Hy  +  y  (1  -  1/Pr)  uuyj  (62) 

Clutter  and  Smith  use  the  following  transformations: 

n  =  (uc/ p*u*x) J  p  dy  (63a) 

and 

x  =  x  (63b) 

They  also  defined  a  stream  function  ty  by 

pru  =  (i|*r)y  and  pvr  =  -('l»r)x  (64) 

The  dimensionless  stream  function  f  is  then  defined  as 
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3f/3n  =  £’  =  u/Uq  (65) 

Then  they  state  the  following  relation  between  f  and  i|>: 

ip  =  (p*n*xue)1/2  f  (66) 

which  is  incorrect.  The  correct  expression  should  read 

<pr  =  (p*y*uex)*/2  /  rf'dn  (67) 

If  r  =  r(x)  only,  then  Equation  67  would  reduce  to  Equation  66  as  given 
by  Clutter  and  Smith;  however,  such  is  not  the  case  since  r(x,y)  = 
rw(x)  +  y  cos  “•  This  error  was  later  recognized  by  Clutter  and  Smith 
[23]  but  not  corrected  until  the  recent  work  of  Jaffe,  Lind,  and  Smith 
[22],  Following  Clutter  and  Smith,  Equation  66  is  assumed,  and  trans¬ 
formation  of  Equation  61  gives 

(C*r)_1  (Crf") 1  +  8(pe/p  -  f’2) 

+  [ (B+l)/2  +  R]ff”  -  x(f'f£  -  f"fx)  =  0  (68) 

where 


and 


C  =  py/peye 
C*  =  p*y*/peue 
B  =  (x/ue)due/dx 

R  =  (x/r)  3r/3x 


(69a) 

(69b) 

(69c) 


(69d) 
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The  no-slip  boundary  condition  is 

£•(0)  =  0  (70) 

Neglecting  mass  transfer  at  the  wall  gives 

fw  =  f(0)  =  0  (71) 

The  remaining  boundary  condition  is 

f»  (n  -*■  «)  =  f'  (•)  »  1  (72) 

if  first-order  vorticity  is  neglected. 

After  defining  the  enthalpy  ratio  g  *  H/H0,  the  energy  equation 
is  transformed  to  thf'  x,n-plane  and  becomes 

(C*r)-1{r[C/Pr)g*  +  (ue2/lle)C(l-l/Pr)f  *f"]  }• 

«  -  [ (B+l)/2  +  R]fg*  +  x(f»gx  -  g*fx)  (73) 

with  the  boundary  conditions 

g(«)  =  1  (74a) 

and  either  g(0)  *  gw(x)  or  g'(0)  ■  g^(x)  (74b) 

Clutter  and  Smith  introduced  the  further  transformations 

<j>'  =  ft  -  1  and  ip  =  g  -  1  (75) 

to  improve  the  accuracy  of  numerical  calculations.  They  also  defined  a 
transverse  curvature  Tcg  by 
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(l/r)3r/3n  =  Tcs  pe/p  (76) 

In  terms  of  these  new  variables,  the  momentum  equation  becomes 
(Cr)'  -  -  Tcs(pe/p)C*"  +  C*0[*'2  +  2<p‘  +  1] 

-  C*N(<j>+n)4>"  +  C*x[(*'+l)<fr'x  -  $»4>x] 

+  (C*(pe/p)x(p*/pe) (u*/ue)2  [ (Y- 1)/y] 

*  [He/u*2  -  1/2]  d(p/p*)/dx}  (77) 

and  the  energy  equation  is 

{(C/Pr)*'  +  (ue2/Ile)  C[ (1-1/Pr)  ($,+l)$"} ' 

-  "  Tcs(pe/e)  UC/Pr)*'  +  (ue2/He)  C[(l-1/Pr)  (♦'♦l)*") 

-  C*N(*+n)*'  +  C*x[(<J»'+l)i|;x  -  *'*x]  (78) 

where  N  »  (0+1) /2  +  R. 

II.  EXTERNAL  VORTICITY  EFFECTS 

Inclusion  of  vorticity  does  not  alter  the  basic  equations  but  only 
the  boundary  conditions.  Hayes  and  Probstein  [12,  page  370]  define  a 
vorticity  index 

0  = 

;e  n*l 

‘wT  f  9iy  (79) 


vorticity  at  the  outer  edge  of  the  boundary  layer 
average  vorticity  across  the  boundary  layer 
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where  ce  is  the  inviscid  vorticity  at  the  outer  edge  of  the  boundary 
layer.  It  should  be  noted  that  the  vorticity  index  defined  by  Equation 
79  will  include  the  effect  of  displacement  when  coupled  with  vorticity 
since  the  properties  are  then  to  be  evaluated  at  effective  body  sur¬ 
face. 

From  Crocco's  vorticity  law  in  terms  of  Van  Dyke's  variables, 
the  inviscid  vorticity  becomes 

Cw  3  (3Ui*/ay*)w  -•  (U*/r*)3U|/3y 

-  -  (U*/r*)  [rwjRTS1’]e  (80) 

where  the  ontropy  gradient  is  obtained  from  Van  Dyke  [14]  as 

-1.6 CM*2  -  l)2  B2sh 

sl(0)  - - 2 - ^“T 

1  (2.8  M2  -  0.4) (2  +  0.4  M„) 

assuming  y  ■  1.4  with  Bsh  =  rn/rs  (see  Equations  20  and  22).  A  second- 
order  expansion  of  the  velocity  leads  to 

u(x,y;e)  =  Uj(x.y)  +  eu2(x,y) 

*  Uj  (x,0)  +  cU2(s,0)  -  ylr^RjTjSlly.Q  (81) 

where  use  was  made  of  Equations  22a,  50a,  and  52a. 

The  vorticity  index  in  the  Clutter-Smith  transformed  plane 

becomes 


Peue 


n=l 

f  p  dy 
o 


Je_ 

Peue 


p»u*  xl 

ue  j 


1/2 


(82) 
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Substituting  liquation  80  into  Equation  82  yields 


fl  »  - 


(«i/V> 

p*  U* 
00  00 


V  R  T  S,'  ■ 

p*u*  x 

■  1/2 

L  ™  Je 

L  ue 

. 

(83) 


Dividing  Equation  81  by  Ug/U^  a  U^(x,0)  eU2(x,0),  substituting 

Equations  63a  and  83,  and  evaluating  the  result  as  y  ■+•  »  yields 


-fM.Jtaw  7  £<>n 

G 


(84) 


for  the  outer  boundary  condition  on  the  momentun  equation  including 
external  vorticity. 

It  should  be  noted  from  Equation  84  that  the  boundary  condition 
f ' («)  must  be  obtained  by  iteration  since  the  integral  is  unknown 
a  priori.  The  numerical  procedure  used  for  solution  of  the  first- 
order  boundary-layer  equations  is  described  in  Appendix  A.  The  pro¬ 
cedure  used  by  Clutter  and  Smith  is  convenient  when  iteration  is 
required,  since  the  basic  method  is  iterative  somewhat  unlike  the 
implicit  finite-diffej snee  scheme  of  Flugge-Lotz  and  Blottner. 

The  vorticity  index  D  can  be  written  as 


«  *  Krwj[(He/u*2)  -  I  (ue/u*)2l  x1/2  (u*/ue)3/2 


(85a) 


where 


K  =  -  (P*/pJ(u*/uJS1'/(Re*)1/2  (85b) 

independent  of  x. 

Thus  one  finds  that  8  'v  x3/2  as  x  -*■  «  for  spherically  blunted 
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cones.  This  index  is  not  to  be  confused  with 

=  rwJ  RjTjSJ  =  [y/Cy-I^VPjS^ 

defined  by  Van  Dyke  [14]  for  the  inviscid  vorticity.  Davis  and 
Flugge-Lotz  [13]  noted  that  flj  x  for  large  x  since.  Sj(0)  = 
constant  and  Pj  -*■  constant  far  downstream  on  most  analytic  bodies  and 
all  spherically  blunted  cones. 

III.  FIRST-ORDER  SLIP  AND  TEMPERATURE  JUMP  EFFECTS 

From  second-order  boundary- layer  theory,  Van  Dyke  [14]  showed 
that  the  si ip- flow  boundary  conditions  are 

u(0]  =  [p/p  (hRT/2) 1/2  3u/3y]w 

and 

T(0)  =  Tw  ♦  (15/8) [u/p  (irRT/2)1/2  3T/3y]w 

when  the  hard-sphere  model  for  the  gas  was  used  to  evaluate  the 
coefficients 

Transforming  the  equations  to  the  x,n-plane  and  expressing  in 
terms  of  convenient  variables  for  inclusion  in  the  Clutter  and  Smith 
numerical  scheme  gives 


f'(n=0) 


PU 


~  Fs  f"(n=0) 


(86a) 


iThe  second  term  in  the  wall  slip  velocity  boundary  conditions 
(Equation  42)  is  of  third  order  (i.e.  of  order  e^)  and  was  not 
considered.  However,  for  the  constant  wall  temperature  conditions 
treated,  the  term  would  have  been  zero  in  any  event. 
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where 


h(n»0)  =  hw 


15  pw  p*_  _  3h 
8  p*u*  p  s  “5? 


(86b) 


F  =  [tt  JL.  u*2  Href  JL  l]  V2  (86c) 

s  \  y-1  2href  hw  u*  Re*  x ) 

Subscript  w  denotes  properties  evaluated  at  the  no-slip  wall  con¬ 
ditions'.  The  fluid  properties  should  be  evaluated  at  the  slip  con¬ 
ditions  as- described  by  Patterson  [24];  however,  to  make  direct  com¬ 
parison  with  the  second-order  results,  the  properties  were  consistently 
evaluated  at  the  no-slip  wall  conditions.  In  recent  calculations  for 
flow  over  a  sharp  cone  at  low  Reynolds  number  conditions,  Mayne, 

Gilley,  and  Lewis  [25],  using  a  modified  version  of  the  Jaffe,  Lind, 
and  Smith  [22]  theory  and  slip-flow  equations  similar  to  those  used 
herein,  evaluated  the  fluid  properties  at  the  slip  conditions.  Un¬ 
published  comparison  of  results  using  fluid  properties  evaluated  at 
both  the  slip  and  no-slip  conditions  showed  small  differences  in  the 
wall  shear  stress. 

The  STJ  boundary  conditions  in  Equations  86  must  be  obtained  by 
iteration  since  the  solution  yields  f"(nB0)  and  3h/3y  /  g' (n=0)  with 
STJ.  As  in  the  case  with  external  vorticity,  inclusion  in  the  Clutter 
and  Smith  numerical  scheme  presents  no  difficulty,  and  the  method  of 
inclusion  is  described  in  Appendix  A. 


IV.  DISPLACEMENT  EFFECT 


In  a  recent  paper,  Van  Dyke  [26]  discussed  the  role  of  the 
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displacement  effect  in  higher-order  boundary- layer  theory.  He  indi¬ 
cated  that  the  first  discussion  of  higher-order  effects  was  given  by 
Prandtl  [27,  page  90]  when  he  described  the  displacement  effects  of  the 
flat-plate  boundary  layer  on  the  uniform  outer  flow  as  follows: 

The  displacement  of  the  stream-lines  by  the  amount  6*  pro¬ 
duces  a  slight  alteration  in  the  potential  flow  which  was  made 
the  basis  of  the  calculations.  Instead  of  the  simple  parallel 
flow,  the  flow  around  a  parabolic  cylinder  of  thickness  26* 
should  be  introduced,  which  would  slightly  alter  the  pressure 
distribution.  The  above  calculation  would  have  to  be  repeated 
for  this  new  pressure  distribution  and  if  necessary  the  process 
repeated  on  the  basis  of  the  new  measure  of  displacement  so 
obtained.  Such  calculations  have  so  far  not  been  performed; 
they  would,  in  any  case,  make  little  difference  in  the  regions 
where  the  calculations  are  usually  applied  in  practice.  They 
would  however,  become  necessary  if  the  transition  to  smaller 
Reynolds  numbers  U^p^L/y^  were  attempted. 

Van  Dyke  [26] ,  in  describing  the  difficulty  in  treating  the  dis¬ 
placement  effect,  noted  the  global  rather  than  the  local  nature  of  the 
effect.  He  termed  the  effect  due  to  classical  boundary- layer  theory 
plus  the  flow  due  to  displacement  as  "one-and-a-half-order- theory." 
However,  the  term  first-order  displacement  effect  will  be  used  herein 
to  describe  the  effect  since  the  choice  of  title  is  arbitrary  and  not 
derived  by  analysis. 

No  modification  to  the  basic  boundary- layer  equations  is 
required  to  treat  the  first-order  displacement  effect.  The  treatment 
of  first-  and  second-order  displacement  effects  is  given  in  Appendix  B, 
whore  the  method  of  treating  the  perturbed  inviscid  outer  flow  is 
described. 
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CHAPTER  IV 

NUMERICAL  RESULTS  AND  DISCUSSION 

Comparisons  of  first-  and  second-order  treatments  of  transverse 
curvature,  vorticity,  slip  and  temperature  jump  (as  a  single  effect), 
and  displacement  will  be  presented  in  this  chapter.  A  sensitive 
quantity  for  such  comparisons  is  the  effect  on  zero-lift  drag,  and  thus 
this  will  be  used  primarily  for  comparison  of-  numerical  results  in  this 
chapter  and  with  experimental  data  in  the  next  chapter.  The  effects  of 
viscosity  law  are  also  easily  shown  by  comparison  of  second-order 
solutions  in  physical  variables  using  a  square-root  viscosity  law  and 
solutions  in  transformed  Levy-Lees  variables  using  the  Sutherland  vis¬ 
cosity  law.1  Unless  otherwise  noted,  all  solutions  used  perfect  gas 
(y  -  1.4),  constant  Prandtl  number  (Pr  ■  0.71),  uniform  wall  tempera¬ 
ture  (Tw/T0  =  0.2  at  Mo#  ■  £  and  0.066  at  M»  ■  18),  Sutherland  viscosity 
law,  and  identical  pressure  distributions  from  an  ideal-gas  blunt  body 
and  characteristics  solution  at  ■  9  or  18. 

For  computation  purposes,  three  sets  of  conditions  or  cases  were 
defined  at  M,,,  =  18  and  five  cases  at  M*  =  9.  The  free-stream  and 

iunpublished  comparisons  between  solutions  in  physical  and  trans¬ 
formed  variables  have  shown  identical  results  for  all  second-order 
effects  except  displacement  and  vorticity.  Therefore,  unless  otherwise 
noted,  the  differences  shown  between  second-order  solutions  were  due 
entirely  to  the  choice  of  viscosity  law. 
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stagnation  conditions  used  in  the  first-  and  second-order  boundary- 
layer  calculations  are  given  in  Tables  I  and  II. 

The  second-order  solutions  for  each  vorticity  and  displacement 
based  upon  displacement  pressure  treatment  were  found  to  be  numerically 
nonunique.  The  solutions  were  affected  by  both  Ay  or  An  step  sizes  and 
n*  or  its  equivalent  in  physical  variables.  It  was  found,  however,  that 
the  sum  of  both  displacement  and  vorticity  effects  was  unique  and  was 
equal  to  the  displacement  speed  treatment  of  the  sum  of  both  effects. 
However,  since  the  displacement  speed  treatment  greatly  magnifies  both 
displacement  and  vorticity  effects  and  since  the  effects  of  An  step  size 
and  n„  were  certainly  significant  but  not  large,  the  second-order 
solutions  shown  were  based  upon  the  displacement  pressure  treatment. 
Moreover,  since  the  sum  of  the  two  effects  is  unique,  no  loss  in 
accuracy  is  involved  when  comparisons  are  made  with  their  combined 
effects. 


I.  TRANSVERSE  CURVATURE  EFFECTS 

A  comparison  of  the  results  of  the  first-  and  second-order 
treatments  of  the  transverse-curvature-induced  friction  drag  at 
M,,  «  9  and  18  is  shown  in  Figures  11  and  12.  At  both  Mach  number  con¬ 
ditions,  the  differences  are  small  between  second-order  solutions  with 
different  viscosity  laws;  however,  the  approximate  treatment  of  Clutter 
and  Smith  predicts  an  increment  about  30  per  cent  less.  This  is  a 
significant  error,  and  the  Clutter  and  Smith  treatment  cannot  be 
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TABLE  I 

EXPANSION  PARAMETER  AND  REYNOLDS  NUMBER 


Case 

M 

00 

e 

Re«, 

Re* 

187 

18.0 

0.0776 

4844.9 

429.1 

185 

18.0 

0.2330 

538.3 

47.7 

182 

18.0 

0.3516 

236.2 

20.9 

95 

9.0 

0.1000 

925.49 

250. 35 

93 

9.0 

0.2401 

160.26 

43.35 

90 

9.0 

0 . 3336 

83.18 

22.5 

91 

9.0 

0.4398 

47.87 

12.95 

92 

9.0 

0.5337 

32.50 

8.79 
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TABLE.  II 

CONDITIONS  USED  IN  THE  BOUNDARY-LAYER  CALCULATIONS 


18.0 

9.0 

He,  ft2/sec2 

4.988xl07 

3.261x10 

lbf/ft* 

146.88 

13.036 

T’.  *R 

8289. 

5263. 

u*2/2href 

0.03923 

0.02565 

Ugai  £t/sec 

10074.8 

7840. 

P„,  lbf/ft2 

0.35172 

0.01244 

p  ,  lbf-sec2/ft* 

1.624X10*6 

2.367x10 

To,.  °R 

126.0 

306.0 

T„,  *R 

547.0 

1053. 

gW 

0.066 

0.20 

Bsh 

0.770 

0.769 
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recommended  for  calculations  involving  transverse  curvature.2 

II.  VORTICITY  INTERACTION  EFFECTS 

Figure  13  shows  the  local  skin- friction  coefficient  predicted 
by  first-  and  second-order  methods.  Good  agreement  exists  for  about 
8.5  nose  radii  from  the  stagnation  point.  Beyond  that  point  the 
first-order  prediction  increases  whereas  the  second-order  result 
remains  approximately  constant  to  the  base  of  the  cone  (x  =  16.4). 

The  effects  of  both  Ax  and  in  step  sizes  were  investigated  for 
the  first-order  solutions.  It  was  found  that  the  solution  was 
sensitive  to  An  but  insensitive  to  Ax  step  sizes.  This  is  explained 
as  follows:  First,  the  transformations  used  by  Clutter  and  Smith  are 
such  that  without  vorticity  effects  the  classical  boundary  layer  was 
substantially  reduced  from  =  6  at  x  »  0  to  n„  =  2  over  much  of  the 
conical  afterbody  (for  longer  bodies  or  bodies  with  more  favorable 
pressure  gradient,  values  of  nB  *  1  were  not  uncommon).  The  variation 
in  'boundary-layer  thickness  was  less  in  the  Levy-Lees  transformed 
plane  than  in  either  the  physical  or  Clutter  and  Smith  transformed 
plane.  Therefore,  for  classical  boundary- layer  solutions,  the  Clutter 
and  Smith  transformation  is  not  optimum  for  bodies  with  strong  favor¬ 
able  pressure  gradients.  Secondly,  and  more  importantly,  the  extent 

Unpublished  comparisons  among  numerical  calculations  of  trans¬ 
verse  curvature  effects  using  a  finite  difference  solution^in  Levy-Lees 
variables,  the  second-order  results  using  the  Davis  and  Fliigge-Lotz 
method,  and  the  first-order  treatment  of  Jaffe,  Lind,  and  Smith  [22] 
showed  that  all  three  methods  gave  essentially  the  same  results  when 
identical  conditions  were  used. 
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Fig.  13  Effects  of  External  Vorticity  on  Skin  Friction  at  M,*,  =-  18,  Tw/T„  =  0.066, 
and  r/c  0.352 


AEDC-TR.  68-148 


AEDC-TR-6B.  148 


of  the  boundary  layer  with  external  vorticity  is  limited  by  a  line  of 
maximum  velocity  »  l/Cl-hg/l^)1/2  (see  Appendix  A).  Along  this 
line  ■  0  and  p^  -*■  *.  This  substantially  reduces  n,,,  below  the  allowed 
classical  solution  value.  It  was  therefore  necessary  occasionally  to 
reduce  the  An  step  size  as  the  solution  proceeded  along  the  body.  The 
data  normal  to  the  surface  were  interpolated  so  that  a  minimum  of  200  n 
points  were  used  to  integrate  the  momentum  and  energy  equations.  The 
solutions  were  accepted  when  similar  calculations  were  made  with  reduced 
step  sizes  and  no  significant  changes  were  observed  in  f"  and  g^. 

Similar  reduction  of  Ax-step  size  significantly  increased  machine  time 
without  affecting  the  solution.  Examination  of  the  f' (n)  profiles 
showed  they  asymptotically  approached  a  linear  variation  as  n  n,,,. 

The  first-order  solutions  shown  with  vorticity  were  examined 
for  step  size  effects  and  are  believed  to  be  unaffected  of  both  Ax  and 
An  step  sizes.  The  divergence  of  solutions  for  x  >  8.5  is  believed  to 
be  a  result  of  the  first-  and  second-order  boundary- layer  treatments 
and  not  an  accuracy  or  step  size  effect.  This  effect  was  observed  for 
all  solutions,  but  the  divergence  in  solutions  increased  as  the  Reynolds 
number  decreased  or  c  increased. 

The  good  agreement  between  the  values  of  Cf  from  the  first-  and 
second-order  treatments  of  the  effects  of  vorticity  for  x  <  8.5  is 
gratifying  and  might  have  been  expected;  however,  to  the  author’s 
knowledge,  this  is  the  first  time  the  essential  agreement  between 
first-  and  second-order  solutions  has  been  demonstrated  for  conditions 
far  from  the  stagnation  region. 
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The  importance  of  the  difference  in  first-  and  second-order 

shear  stress  is  shown  in  Figures  14  and  15.  At  M  *  9  the  first-  and 

00 

second-order  solutions  were  in  excellent  agreement  for  e  <  0.2  and 
departed  rapidly  as  e  increased.  For  e  =  0.4  the  first-order  predic¬ 
tion  was  approximately  twice  the  second-order  increment.  The  large 
difference  in  this  the  largest  higher-order  drag  component  is  dis¬ 
appointing,  and  moreover,  it  will  be  shown  by  comparison  with  experi¬ 
mental  data  that  both  methods  overpredict  the  total  zero- lift  drag  for 
e  >  0.2.  The  trends  of  the  results  at  =  18  were  similar  to  those 
at  =  9,  but  the  solutions  tended  to  separate  earlier  and  one  of  the 
first-order  solutions  fell  below  the  second-order  solution  curve.  The 
range  of  e  was  chosen  for  numerical  solution  to  cover  generally  the 
range  of  available  experimental  data,  and  thus  the  ranges  differed  for 
the  two  Mach  numbers. 

Figure  16  shows  a  convenient  quantity  introduced  by  Davis  and 

•  •  • 

Flugge-Lotz  for  a  stretched  nondimens ionali zed  displacement  thickness. 

The  quantity  (which  they  denoted  as  DISP)  is  defined  here  as  5g  =  5*/e 

and  thus  is  independent  of  the  expansion  parameter.  Comparisons  of  the 

results  shown  on  Figure  16  indicate  the  effects  of  viscosity  law  on 

the  two-dimensional  displacement  thickness,  the  effects  of  transverse 

curvature,  and  the  effects  of  first-order  vorticity  interaction  on  the 

displacement  thickness.  Using  a  square-root  viscosity  law  gave  a 

slightly  larger  6*  than  the  use  of  Sutherland's  law  predicted.  The 

axisymmetric  5*  including  TVC  was  slightly  below  the  two-dimensional 

result.  A  dramatic  effect  occurred  when  the  expression  for  the 
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acd, 


Fig.  14 


First-  and  Second-Order  Vorticity-lndoced  Friction  Drag  Predictions  at 
=  9,  Tw/T„  =  0.2,  and  r„/rb  =  0,3 
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Formula 


Method 


|i  Law 

- 2- Dimensional  Square  Root  Davis  Physical  Variables 

- 2-Dimensional  Sutherland  Davis  and  Levy- Lees  Variables 

-  Axisymmetric  +  TVC  Sutherland  Clutter  and  Smith 

- Axisymmetric  +  Vorticity  Sutherland  Present  Results 
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displacement  thickness  6*.  which  is  derived  in  Appendix  C  and  includes 
vorticity  interaction,  was  compared  with  the  classical  forms.  Includ- 
ing  vorticity  produced  a  negative  displacement  thickness  over  the  nose 
(whereas  the  values  calculated  using  classical  theory  were  slightly 
positive),  substantially  increased  the  rate  of  growth  in  the  region  of 
strong  favorable  pressure  gradient,  and  caused  a  rapid  decline  in  dis- 
placement  thickness  over  much  of  the  conical  afterbody  where  the 
effects  of  vorticity  on  the  velocity  profile  and  wall  shear  were  most 
important.  The  behavior  of  6*  indicates  the  difficulty  of  realistically 
and  consistently  treating  the  combined  effects  of  displacement  and 
vorticity  interactions  by  a  first-order  method  without  simultaneously 
treating  the  interaction  of  the  entropy  and  boundary  layers.  Of  course, 
in  the  second-order  theory  this  interaction  was  not  considered  and  its 
region  of  applicability  was  also  limited.  These  limitations  can  also 
be  seen  in  the  next  chapter  from  comparisons  of  the  numerical  results 
with  some  experimental  data. 

III.  DISPLACEMENT  EFFECTS 

The  first-order  displacement- induced  pressure  distributions  over 
the  spherically  blunted  cone  are  shown  in  Figures  17  and  18.  The  large 
effect  of  displacement  on  the  pressure  distributions  for  the  ranges  in 
Reynolds  number  considered  is  evident  from  the  figures.  The  bump  in 
the  pressure  distribution  for  the  lowest  Reynolds  number  or  largest  e 
considered  at  M*  *  9  can  be  seen  in  Figure  18.  The  bump  is  caused  by 
the  rapid  increase  in  displacement  thickness  in  the  region  x  «  1.0  to 
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2.0.  As  described  in  Appendix  B,  the  curve  fit  of  the  effective  body 
in  the  region  of  rapid  increase  in  displacement  thickness  was  poorer 
since  inflection  points  were  not  allowed  in  the  effective  body  curve 
fits  over  this  region.  However,  allowing  inflection  points  in  the 
curve  fits  would  only  have  caused  larger  bumps  and  induced  secondary 
shocks  which  would  have  further  altered  the  pressure  distributions 
over  the  conical  afterbody. 

First-  and  second-order  displacement-induced  pressure  and 
friction-drag  increments  are  shown  in  Figures  19  and  20.  Perturbed 
inviscid  outer  flow  solutions  were  obtained  for  four  values  of  e 
at  Mw  =  9  and  three  values  at  M,,,  ■  18.  The  two  predictions  for  the 
displacement-induced  pressure  drag  differed  only  by  the  effect  of  the 
pressure  gradient  over  the  nose,  where  the  effect  was  included  in  the 
second-order  but  not  in  the  first-order  treatment. 

The  differences  between  the  two  treatments  of  displacement- 
induced  friction  drag  were  larger,  and  the  results  are  more  interest¬ 
ing.  In  the  first-order  treatment,  only  the  effect  of  displacement 
of  the  inviscid  outer  flow-field  pressure  was  included.  In  the  second- 
order  solutions,  however,  the  effects  of  inviscid  vorticity  on  the 
outer  edge  velocity  and  temperature  were  also  present. 

Caution  must  be  used  when  trying  to  generalize  the  results 
shown  in  Figures  19  and  20.  First,  the  second-order  results  for 
displacement-induced  pressure  and  friction-drag  increments  were  not 
simply  linearly  and  quadratically  dependent,  respectively,  on  e. 

(It  is  interesting  to  note,  however,  that  the  first-order  friction-drag 


69 


AEDC-TR-68-148 


AEDC-T  R •68-146 


0.16 


Order 

Variables 

Li  Law 

1st 

Clutter  and  Smith 

Sutherland 

2nd 

Levy-Lees 

Sutherland 

2nd 

Physical 

Square  Root 

Fig.  20  First*  and  Second-Order  Displacement-Induced  Pressure  and  Friction-Drag  Components 
at  Mo,,  =  18,  Tw/T0  =  0.066,  and  rn/rb  =  0.3 
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results  were  quadratically  dependent  on  e.)  Secondly,  for  other  bodies 
under  different  conditions,  negative  second-order  displacement- induced 
friction  drag  has  been  found.  The  effects  of  geometry  (both  surface 
slope  and  length),  wall-to-stagnation  temperature  ratio,  and  viscosity 
law  used  have  the  strongest  influence  on  second-order  wall  shear  stress 
distribution.  For  the  conditions  treated  here,  the  effects  of  wall-to- 
stagnation  temperature  ratio  caused  small  differences  in  the  induced- 
drag  increments  for  both  first-  and  second-order  treatments. 

It  is  emphasized  again  that  the  second-order  displacement-induced 
friction-drag  results  shown  in  Figures  19  and  20  were  not  unique.  The 

results  shown  in  Levy-Lees  variables  were  obtained  with  an  edge  n  =  6. 

00 

Two  conditions  were  recalculated  using  =  9,  and  the  displacement- 
induced  friction  drag  was  increased  about  7  per  cent  without  affecting 
the  displacement -induced  pressure  drag.  Therefore,  the  excellent  agree¬ 
ment  between  the  second-order  solutions  in  physical  and  Levy-Lees 
variables  was  to  some  extent  fortuitous  .since  neither  solution  was 
independent  of  the  outer  matching  point  location. 

IY.  SLIP  AND  TEMPERATURE  JUMP  EFFECTS 

Comparisons  of  first-  and  second-order  wall  slip  velocity  and 
temperature  jump  for  one  of  the  conditions  at  ■  9  are  shown  in 
Figures  21  and  22.  The  first-order  slip  velocity  and  temperature  was 
lower  than  the  second-order  results  by  about  15  to  20  per  cent.  Figure 
23  shows  a  comparison  of  the  first-  and  second-order  local  (total) 
skin-friction  coefficient.  The  differences  in  slip  velocity  and 
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Fig.  21  First-  and  Second-Order  Slip  Velocity  at  M«,  =  9,  Tw/T0  =  0.2,  and  e  =  0.24 
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Fig.  23  Effects  of  Slip  and  Temperature  Jump  an  Skin  Friction  at  M^,  =  9,  Tw/T0  =  0.2, 
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temperature  caused  a  negligible  effect  on  the  wall  shear  stress  dis¬ 
tribution.  Moreover,  for  the  lowest  Reynolds  number  condition  con¬ 
sidered  (Re in.  <  400),  the  separate  effects  of  slip  and  temperature 
jump  were  indeed  small  in  comparison  with  other  effects. 

The  slip-and-temperature-jump-induced  friction-drag  results 
are  shown  in  Figures  24  and  25,  At  the  highly  cooled  wall  conditions 
at  =  18,  the  second-order  solution  in  physical  variables  was  step 
size  dependent,  and  an  acceptable  solution  was  not  obtained  for  this 
effect  at  these  conditions  using  physical  variables.  In  either  the 
Clutter  and  Smith  or  Levy-Lees  transformed  planes,  the  highly  cooled 
wall  presented  no  difficulty,  and  stable  solutions  were  obtained  with 
■reasonable  (An  ■  0.025  or  0.05)  step  sizes.  Both  first-  and  second- 
order  treatments  evaluated  the  properties  in  the  slip  velocity  and 
temperature  at  the  no-slip  wall  conditions.  The  differences  between 
first-  and  second-order  treatments  were  larger  for  this  separate  effect. 
However,  since  this  drag  component  was  small  compared  with  other  second- 
order  effects,  the  difference  between  the  two  treatments  was  not  an 
important  influence  on  the  total  drag  prediction  for  e  <  0.2. 

A  comparison  of  the  first-  and  second-order  incremental  friction- 
drag  components  at  M,„  =  18  is  shown  in  Figure  26.  The  nonlinearities 
of  the  first-order  vorticity-induced  friction  drag  and  the  second- 
order  displacement-induced  friction-drag  components  arc  clearly  visible 
as  well  as  the  empirically  linear  variation  of  the  first-order  dis¬ 
placement  and  TVC  effects.  The  second-order  results  in  physical  vari¬ 
ables  using  a  square-root  viscosity  law  are  not  shown  since  they  were 
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Fig.  24  First-  and  Second-Order,  Slip-and-Temperature-Jump-Induced  Friction  Drag  atMK  =  9, 
Tw/Te  =  0-2,  and  rn/rb  =  0.3 
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Order 


Variables 


M  Law 

- 2nd  Levy-Lees  Sutherland 
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Fig.  25  First-  and  Second-Order,  SI  ip-and-Temperature-Jump-Induced  Friction  Drag  at  M«,  =  18, 
Tw/Te  =  0.066,  and  rn/rb  =  0.3 
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in  good  agreement  with  the  second-order  results  presented.  Displayed 
in  this  way,  it  is  obvious  that  for  the  conditions  considered  all  of 
the  separate  higher-order  drag  components  can  be  classified  as  either 
major  (TVC,  displacement  and  vorticity)  or  minor  (STJ  and  LC)  effects. 
However,  when  the  effects  are  coupled  in  the  first-order  treatment,  no 
simple  classification  exists.  The  results  of  such  a  comparison  at 

*  9  showed  qualitative  agreement  with  the  results  at  M,,,  =  18;  there¬ 
fore,  the  *  9  data  are  not  presented. 

V.  COUPLED  EFFECTS 

The  first-order  treatment  of  higher-order  effects  on  the  skin- 
friction  coefficient  at  M„  «  9  is  shown  in  Figure  27.  The  results  are 
shown  for  a  large  value  of  the  expansion  parameter  e  =  0.533  in  order 
to  magnify  some  of  the  effects  and  clearly  establish  the  strong  coupling 
influence  of  small  individual  effects. 

The  individual  effects  of  transverse  curvature  (TVC) ,  vorticity 
(vort),  displacement  (disp) ,  and  slip  and  temperature  jump  (STJ)  on  the 
skin-friction  coefficient  are  seen  by  comparison  with  the  classical 
first-order  axisymmetric  result.  The  bumps  in  the  curves  including 
displacement  were  caused  by  a  displacement-induced  compression  region 
over  the  effective  body  and  the  resulting  displacement  of  the  inviscid 
outer  flow.  Although  at  the  conditions  considered  in  Figure  27,  each 
separate  higher-order  effect  had  a  significant  influence  on  the  skin- 
friction  distribution,  the  separate  effect  of  vorticity  was  clearly 
dominant.  The  strong  effect  of  coupling  displacement  and  vorticity  can 
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Fig.  27  Separate  and  Coupled  Higher-Order  Effects  an  Skin-Frictian  Coefficient,  =  9, 
Friction  Coefficient  at  =  9.  T^/T_  =  0.2.  and  t  =  0.533 
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be  seen  by  comparison  of  the  results  of  these  combined  effects  with  the 
sum  of  the  independent  effects.  Coupling  displacement  with  vorticity 
approximately  doubled  the  vorticity-induced  skin-friction  coefficient. 
The  displacement-induced  pressure  strongly  affected  the  vorticity  index 
fl  which  in  turn  controlled  the  vorticity  effect.  Coupling  TVC  with 
displacement  and  vorticity  increased  the  coupled  effects  of  the  latter 
two  by  about  the  increment  due  to  transverse-curvature-induced  skin 
friction  only.  The  most  surprising  result  was  found  when  slip  and 
temperature  jump  effects  were  coupled  with  displacement,  vorticity,  and 
transverse  curvature.  Near  the  end  of  the  body,  the  effects  of  STJ 
were  larger  than  the  combined  effects  of  displacement  and  transverse 
curvature!  The  trends  clearly  indicate  that  if  the  body  were  longer 
the  effects  of  slip  and  temperature  jump  would  offset  not  only  the 
effects  of  displacement  and  transverse  curvature  but  vorticity  as  well. 
In  Chapter  III  it  was  shown  that  f 1 (n®0)  «  f'Ch^O).  Vorticity  increased 
f"  which  in  turn  increased  the  slip  velocity  f ' (n~0)  which  reduced  the 
velocity  gradient  r'(nB0).  This  coupling  between  vorticity  and  STJ  can 
be  seen  by  comparison  of  the  almost  linear  increase  in  the  vorticity- 
induced  skin-friction  coefficient  with  the  almost  linear  decrease  in 
the  combined  total  effects  including  slip. 

The  effects  of  coupling  first-order  displacement  and  vorticity 
effects  can  be  seen  in  Figures  28  and  29.  The  vorticity  index  Sl(x) 
is  shown  for  the  conditions  at  M,,  =  9  and  18  both  with  and  without  the 
displacement  effect.  As  the  Reynolds  number  decreases,  e  increases  for 
both  Mach  numbers.  The  dependence  of  8(x)  on  the  pressure  can  be  seen 
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as  follows:  In  Chapter  III  it  was  shown  that 

fl(x)  *  [  (Hg/u*2)  -  I  Cue/uJ2]  (u*/ue)3/z 

when  x  is  held  fixed  and  the  pressure  allowed  to  vary.  From  the 
perfect  gas  law  and  energy  equation 

p  =  PRT  »  III  ph  =  Ilip (He  -  ue2/2) 

Along  a  streamline  p  «  pY  ,  and  also 


Therefore,  one  finds  that 


where  the  remaining  terms  are  all  positive  and  can  be  neglected  in 
comparison  with  the  three  leading  terms  for  p/p^  <  0.1. 

As  the  Reynolds  number  is  reduced  or  e  increases,  the  displace¬ 
ment  thickness  increases  which  produces  a  larger  displacement- 
induced  pressure  (see  Figures  17  and  1-8,  pages  67  and  68).  Thus,  when 
the  displacement  and  vorticity  effects  are  coupled,  the  displacement- 
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induced  pressure  increases  the  vorticity  effect  by  increasing  D(x). 

Because  of  the  hotter  wall  at  M«,  =  9  (TW/TQ  =  0,2  at  M*  *  9, 
whereas  Tw/T0  =  0.066  at  =  18),  the  displacement  thickness  was 
positive  over  the  spherical  nose  rather  than  negative  as  in  the  M  = 

18  cases.  The  rapid  increase  in  displacement  thickness  at  *  9  for 
0.8  <  x  <  2.5  caused  bumps  in  the  S2(x)  curves  which  were  not  observed 
at  the  M,,,  =  18  highly  cooled  wall  condition.  It  should  be  noted  that 
e  =  0.33  for  Case  90  with  a  bump  and  e  =  0,35  for  Case  182  without  a 
bump,  and  thus  there  was  a  region  of  overlap  for  comparison. 

Although  the  coupling  between  displacement  and  vorticity  was  not 
as  dramatic  as  the  coupling  between  vorticity  and  STJ,  the  effect  was 
important.  Since  the  first-order  treatment  always  produces  an  increase 
in  displacement- induced  pressure,  the  effects  of  displacement  and 
vorticity  coupling  would  be  expected  to  increase  with  increasing  wall- 
to-stagnation  temperature  ratio  and  increasing  body  length. 

The  results  for  the  coupled  higher-order  effects  were  clearly 
not  linearly  independent  as  were  the  second-order  effects  of  TVC, 
longitudinal  curvature,  vorticity,  and  slip  and  temperature  jump 
according  to  Van  Dyke’s  theory.  The  obvious  result  here  was  that  small 
separate  higher-order  effects  had  a  strong  influence  when  coupled  with 
other  higher-order  effects,  and  a  simple  linear  combination  of  each 
first-order  effect  was  not  justified  and  could  lead  to  erroneous  results, 

VI.  RANGE  OF  APPLICABILITY 

The  range  of  applicability  of  the  first-  and  second-order 
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theories  must  finally  be  determined  by  comparison  with  experimental 
data;  however,  it  is  instructive  to  consider  the  range  of  applicability 
predicted  from  the  numerical  results. 

The  first-order  boundary  layer-to-shock  layer  thickness  ratio 
over  the  cone  at  a  18  is  shown  in  Figure  30.  The  boundary- layer 
thickness  was  defined  to  be  the  distance  normal  to  the  wail  where  f*  = 
u/ue  =  0,1)95  and  y^  is  the  first-order  (nonperturbed)  inviscid  shock- 
layer  thickness  from  the  blunt  body  and  characteristics  solutions.  The 
variation  of  <5/ysh  was  small  over  the  entire  body,  as  shown  in  Figure 
30,  and  a  region  is  indicated  between  e  =  0.11  and  0.16  where  <5/ysh  s 
1.  Beyond  this  region  the  shock  layer  is  fully  viscous  and  boundary- 
layer  theory  is  not  applicable. 

The  values  of  x  shown  in  Figure  30  correspond  to  some  of  the 
stations  along  the  body  used  in  the  first-order  boundary- layer  cal¬ 
culations  in  Clutter  and  Smith  variables.  The  boundary  layer-to-shock 
layer  thickness  ratio  is  minimum  at  the  stagnation  point  (x  =  0),  a 
maximum  at  the  sphere-cone  tangent  point  (x  =  1.409),  and  the  variation 
over  the  entire  body  is  rather  small  (about  30  per  cent). 

Since  the  boundary-layer  thickness  was  arbitrarily  defined,  the 
effect  of  the  definition  of  the  boundary- layer  thickness  on  5/ysh  is 
shown  in  Figure  31.  The  results  of  the  variation  f^  =  0.90  to  0.999  are 
shown  at  the  location  x  =  3.54.  The  results  for  f£  =  0,90  and  0.95  are 
shown  only  for  interest  since  these  values  are  too  low  for  serious  con¬ 
sideration.  The  range  in  c  is  not  significantly  affected  by  the 
variation  f£  =  0.990  to  0.999,  and  thus  the  value  of  0,995  appears  to  be 
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Fig.  30  Boundary  Layer-to-Shock  Layer  Thickness  Ratio  Over  the  Cone  at  M, 
and  Tw/T0  =  0.66 
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Fig.  31  Effects  of  the  Definition  of  Boundary-Layer  Thickness  on  the  Range  of 
Applicability  of  Boundary-Layer  Theory 
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a  reasonable  choice  for  the  definition  of  boundary -layer  thickness. 

VII.  COMPUTING  TIME  REQUIREMENTS 

In  concluding  this  chapter  on  the  discussion  of  numerical  results, 
a  comparison  is  given  of  computing  time  required  for  the  first-  and 
second-order  boundary- layer  calculations.  The  times  required  for  the 
calculations  on  a  Control  Data  Corporation  1604-B  digital  computer  are 
shown  in  Table  III.  The  time  required  for  the  first-order  treatment 
including  the  vorticity  effect  increased  as  the  Reynolds  number  de¬ 
creased  (or  as  c  increased)  for  each  of  the  free-stream  Mach  numbers 
considered.  Therefore,  the  times  shown  in  Table  III  are  valid  only  for 
one  set  of  conditions  (Case  185),  and  the  data  are  shown  primarily  for 
comparison  of  the  machine  time  required  for  the  first-  and  second-order 
boundary-layer  calculations . 3 

From  the  results  shown  in  Table  III,  the  implicit  finite- 
difference  method  in  physical  variables  or  Levy-Lees  variables  was 
approximately  three  times  faster  than  the  first-order  solution  in  Clutter 
and  Smith- variables  for  treating  the  classical  boundary- layer  problem. 

It  should  be  noted  also  that  the  number  of  stations  used  in  the  finite- 
difference  solutions  was  approximately  ten  times  the  number  used  in  the 
first-order  treatment.  Including  the  second-order  effects  in  the  finite- 
difference  solutions  approximately  doubled  the  computing  time  required 

^The  computing  time  required  for  the  blunt  body  and  characteris¬ 
tics  solutions  (approximately  20  minutes  for  each  case  considered)  was 
in  addition  to  the  time  shown  in  Table  III  for  the  boundary- layer 
calculations. 
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TABLE  III 

COMPUTING  TIMES  FOR  FIRST-  AND  SECOND-ORDER 
TREATMENTS  OF  HIGHER-ORDER  BOUNDARY-LAYER  EFFECTS0 


Effect 

First-Order 

Clutter-Smith 

Variables 

Second-Order 

Physical 

Variables 

Second-Order 
Levy- Lees 
Variables 

Computing  time,  minutes*5 

Classical*- 

75 

25 

26 

1YC 

75 

50 

55 

STJ 

75 

50 

55 

disp 

75 

50 

55 

vort 

100 

50 

55 

disp  +  vort 

120 

— 

-- 

disp  ♦  vort  +  TVC 

165 

mi  a* 

-- 

disp+vort+TVC+STJ 

180 

50 

55 

Supplementary  Data 


Stations  along  the 
body 

34 

336 

390 

nm  (maximun) 

3 

— 

6 

An 

0,025 

--- 

0.05 

AN 

ate  mi 

0.02 

aCase  185,  b\m  =  18,  e  =  0.233 

Wimc  required  on  a  Control  Data  Corporation  1604-B  digital 
computer  (computing  speed  approximately  60  per  cent  of  IBM  7094)  with 
an  on-line  printer. 

Calculations  for  axisymmetric  body  without  any  higher-order 
effects. 
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for  the  classical  boundary -layer  problem. 

For  engineering  purposes,  one  would  usually  be  interested  in 
calculations  which  include  all  of  the  significant  effects.  In  both 
first-  and  second-order  treatments,  at  least  two  complete  boundary- 
layer  calculations  were  required  since,  using  the  methods  described 
herein,  the  effects  of  displacement  could  not  be  obtained  with  a  single 
calculation.  Therefore,  for  the  conditions  shown  in  Table  III,  a 
minimum  of  225  minutes  was  required  for  the  first-order  treatment 
including  all  effects  and  75  or  81  minutes  was  required  for  the  second- 
order  treatment  in  physical  or  Levy-Lees  variables,  respectively. 

The  substantial  savings  in  computing  time  using  the  implicit  finite- 
difference  method  makes  this  scheme  attractive  for  other  applications 
(for  example,  the  chemically  reacting  boundary- layer  problem). 


92 


AEDC-TR.  68-148 


CHAPTER  V 

COMPARISON  OF  NUMERICAL  RESULTS  AND  EXPERIMENTAL  DATA 

Results  of  the  first-  and  second-order  treatments  will  be  com¬ 
pared  with  some  available  experimental  data  on  pressure  and  heat- 
transfer  distributions  and  zero-lift  drag  of  a  spherically  blunted  cone 
at  M,,,  ®  9  and  18.  The  experimental  data  were  obtained  from  Whitfield 
and  Griffith  [1]  and  [28],  Lewis  [29],  and  Griffith  and  Lewis  [30]. 

Some  previous  numerical  results  of  Lewis  and  Whitfield  [7]  will  also  be 
compared  with  the  numerical  results  from  the  first-  and  second-order 
treatments  and  the  experimental  data. 

I.  PRESSURE  DISTRIBUTIONS 

The  displacement-induced  pressure  over  a  spherically  blunted  cone 
at  Mm  =  20  is  shown  in  Figure  32.  The  first-order  iterated  results  of 
Lewis  and  Whitfield  [7]  and  the  second-order  results  of  Marchand, 

Lewis,  and  Davis  [18]  at  Mw  ■  18  are  compared  with  the  experimental  data 
of  Lewis  [29]  and  Griffith  and  Lewis  [30],  It  is  obvious  that  the 
trends  of  the  experimental  results  were  not  predicted  by  either  first- 
or  second-order  treatments.  The  nonlinear  character  of  the  first- 
order  results  of  Lewis  and  Whitfield  was  due  to  the  effects  of  iterating 
the  inviscid  outer  flow  and  the  viscous  boundary- layer  flow  fields.  The 
effect  of  this  iteration  was  always  to  reduce  the  displacement- induced 
pressure. 
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Fig.  32  Displacement-Induced  Pressure  Over  a  Spherically  Blunted  Cone  at  Moo 
Tw/T0  =  0.066,  6C  =  9  degrees,  and  rn/rj,  =  0.3 
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The  poor  agreement  between  numerical  results  and  experimental 
pressure  data  is  not  understood.  The  experimental  data  lay  in  a  range 
where  the  second-order  theory  should  have  been  applicable.  Moreover , 
the  experimental  data  and  numerical  results  have  been  checked  extensively, 
and  the  errors  due  to  approximations  in  the  numerical  treatments  and 
experimental  scatter  are  believed  to  be  small  deviations. 

II.  HEAT-TRANSFER  DISTRIBUTIONS 

Figure  33  shows  a  comparison  of  first-  and  second-order  heat 
transfer  to  a  sphere-cone  at  *  18.  In  contrast  with  the  pressure 
data  discussed  above,  the  first-  and  second-order  results  were  in 
excellent  agreement  with  the  experimental  data  of  Griffith  and  Lewis  [30]. 
Because  of  the  limited  experimental  data  available,  caution  must  be  used 
in  interpreting  the  comparison  with  the  numerical  results. 

III.  ZERO-LIFT  DRAG 

The  final  and  most  instructive  comparisons  between  first-  and 
second-order  numerical  results  and  experimental  data  are  given  in  Figures 
34  through  41.  The  total  drag  predicted  by  the  two  treatments  is  com¬ 
pared  with  the  experimental  data  of  Whitfield  and  Griffith  [1]  and  [28], 

The  results  from  the  first-order  treatment  of  the  higher-order 
effects  on  zero-lift  drag  are  given  in  Table  IV.  Comparison  of  the  sum 
of  separate  effects  with  the  results  including  the  coupled  effects  shows 
that  the  results  were  not  linearaly  independent,  but  strong  coupling  was 
involved  especially  at  the  highest  e.  The  relatively  good  agreement 
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Fig.  33  Heat  Transfer  to  a  Spherically  Blunted  Cone  at  M«  =  18,  Tw/T0  =  0.066, 
Be  =  9  degrees,  and  rn/r(,  =  0.3 
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Fig.  34  Predictions  of  the  Drag  of  a  Spherically  Blunted  Cone  at  MK  =  18, 

Tw/To  =  0.066,  and  rn/rt  =  0.3  Using  Second-Order  Boundary-Layer 
Theory  in  Physical  Variables 
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Fig.  36  Predictions  of  the  Drag  of  a  Spherically  Blunted  Cone  at  =  18, 
Tw/T„  =  0.066,  and  rn/rb  =  0.3  Using  First-Order  Boundary-Layer 
Theory  in  Clutter  and  Smith  Variables 
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Numerical  Results 


Van  Dyke  -  Davis  and  Fliigge-Lotz 
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Fig.  37  Comparison  of  First- and  Second-Order  Predictions  of  the  Drag  of  a 
Spherically  Blunted  Cone  at  M,*  =  18,  Tw/T0  =  0.066,  and 
rn/rb  =  0.3 
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Fig.  38  Predictions  of  the  Drag  of  a  Spherically  Blunted  Cone  at  Mm  =  9, 
Tw/T0  =  0*2,  and  r„/rj,  =  0.3  Based  an  Second-Order  Boundary- 
Layer  Theory  in  Physical  Variables 
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Fig.  41  Comparison  of  First-  and  Second-Order  Boundary-Layer  Theory 

Predictions  of  the  Drag  of  a  Spherically  Blunted  Cone  at  =  9, 
Tw/T„  =  0.2,and  rn/rb  =  0.3 
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TABLE  IV 

FIRST-ORDER  RESULTS  OF  HIGHER-ORDER  EFFECTS  ON  ZERO-LIFT  DRAG 


Case 

e 

Effectsa 

C°f 

% 

C°Total 

187 

0.0776 

TVC 

0.0550 

0.0550 

0.1692 

187 

0.0776 

STJ 

0.0514 

0.1142 

0.1656 

187 

0.0776 

vort 

0.0555 

0.1142 

0.1697 

187 

0.0776 

vort+TVC 

0.0584 

0.1142 

0.1726 

187 

0.0776 

disp 

0.0552 

0.1223 

0.1775 

187 

0.0776 

disp+vort 

0.0600 

0.1223 

0.1823 

187 

0.0776 

disp+vort+TVC 

0.0633 

0.1223 

0.1856 

187 

0.0776 

disp+vort+TVC+STJ 

0.0630 

0.1223 

0.1853 

185 

0.2530 

TVC 

0.1851 

0.1142 

0.2993 

185 

0.2530 

STJ 

0.1529 

0.1142 

0.2672 

185 

0.2530 

vort  . 

0.1915 

0.1142 

0.3057 

185 

0.2530 

vort+TVC 

0.2206 

0.1142 

0.3348 

185 

0.2530 

vort+TVC+STJ 

0.2106 

0.1142 

0.3248 

185 

0.2530 

disp 

0.1878 

0.1437 

0.3315 

185 

0.2530 

disp+vort 

0.2558 

0.1437 

0.3995 

185 

0.2530 

disp+vort+TVC 

0.2873 

0.1437 

0.4310 

185 

0.2530 

disp+vort+TVC+STJ 

0.2856 

0.1437 

0.4293 

182 

0.3516 

Classical^5 

0.2335 

0.1142 

0.3477 

182 

0.3516 

TVC 

0.3018 

0.1142 

0.4160 

182 

0.3516 

STJ 

0.2295 

0.1142 

0.3437 

182 

0.3516 

vort 

0.3497- 

0.1142 

0.4639 

182 

0.3516 

vort+TVC 

0.4268 

0.1142 

0,5410 

182 

0.3516 

disp 

0.3118 

0.1667 

0.4785 

182 

0.3516 

disp+vort 

0.5628 

0.1667 

0.7295 

182 

0.3516 

disp+vort+TVC 

0.6304 

0.1667 

0.7971 

182 

0.3516 

disp+vort+TVC+STJ 

0.6229 

0.1667 

0.7896 

95 

0.1000 

TVC 

0.0811 

0.1667 

0,1977 

95 

0.1000 

vort 

0.0819 

0.1667 

0.1985 

95 

0.1000 

disp 

0.0798 

0,1320 

0.2118 

95 

0.1000 

disp+vort 

0.0905 

0.1320 

0.2225 

95 

0.1000 

disp+vort+TVC 

0.0973 

0.1320 

0.2293 

95 

0.1000 

disp+vort +TVC+STJ 

0.0961 

0.1320 

0.2281 

93 

0.2401 

TVC 

0.2185 

0.1166 

0.3351 

93 

0.2401 

STJ 

0. 1724 

0.1166 

0.2890 

93 

0.2401 

vort 

O'.  2464 

0. 1166 

0.3630 
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between  the  first-  and  second-order  results  of  the  total  drag  predictions 
shown  in  Figures  37,  page  100,  and  41,  page  104,  should  be  considered 
fortuitous  in  view  of  the  strong  effects  of  first-order  coupling. 

On  Figures  36,  page  90,  and  40,  page  103,  the  curves  labeled 
C°total  were  obtained  by  summation  of  the  separate  effects  treated 
independently  (that  is,  no  coupling  of  effects  was  considered).  Separate 
curves  labeled  Q)vort+diSp+TVC+STJ  are  shown  which  included  all  effects 
coupled  simultaneously.  For  both  Mach  number  conditions,  the  results 
show  excellent  agreement  between  the  uncoupled  and  coupled  drag  pre¬ 
dictions  for  e  -  0.1  and  somewhat  surprisingly  good  agreement  for  e  >  0.1 
considering  the  strong  coupling  effpets  on  the  skin-friction  coefficient 
shown  in  Figure  27,  page  81. 

For  M0  =  18  and  e  <  0.15,  good  agreement  was  found  between 
* 

numerical  results  and  the  available  experimental  data,  but  for  e  >  0.2  the 
agreement  was  poor.  Recalling  that  in  the  second-order  theory  it  was 
assumed  that  e  <<  1,  one  should  expect  that  as  e  increases  at  some 
value  the  theory  will  no  longer  be  applicable.  For  conditions  considered 
here,  this  point  appears  to  be  near  c  «  0,15.  From  this  comparison  the 
range  of  applicability  of  the  theory  is  consistent  with  the  range  pre¬ 
dicted  from  consideration  of  boundary  layer* to-shock  layer  thickness 
ratio,  <5/ysh. 

The  previous  first-order  results  of  Lewis  and  Whitfield  [7]  are 
also  shown  for  comparison  in  Figures  37  and  41.  The  apparently  good 
agreement  between  their  prediction  and  experimental  data  was  simply 
fortuitous  since  only  approximate  transverse  curvature  (Clutter  and  Smith) 
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and  displacement  (iterated  inviscid-viscous  flow  fields)  were  included „ 

The  comparisons  of  transverse-curvature- induced  friction  drag  showed  that 
the  approximate  Clutter  and  Smith  treatment  led  to  an  error  of  about  30 
per  cent  in  that  component*  Also,  iteration  of  the  inviscid-viscous  flow 
fields  reduced  the  displacement-induced  pressure  and  friction  drag.  In  a 
first-order  sense,  it  can  be  argued  that  iteration  was  allowable  until 
there  was  negligible  change  in  p(x).  The  fact  remains,  however,  that  in 
the  earlier  work  a  sizable  error  existed  in  the  transverse  curvature  term 
and  the  important  contribution  of  external  vorticity  was  not  considered. 

It  is  interesting  to  note  that  the  improved  theoretical  model  and  numerical 
technique  gave  results  in  poorer  agreement  upon  comparison  of  "theoreti- 
cal,,  results  and  experimental  data.  However,  the  application  of  highly 
approximate  schemes  outside  the  range  whore  they  have  been  extensively 
compared  with  experimental  data  can  lead  to  erroneous  results. 
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CHAPTER  VI 

CONCLUSIONS  AND  SUMMARY 

Based  upon  comparison  of  first-  and  second-order  treatments  of 
transverse  curvature,  vorticity,  displacement,  and  slip  and  temperature 
jump,  and  comparisons  of  numerical  results  with  experimental  data,  the 
following  conclusions  are  drawn: 

1.  The  approximate  treatment  of  transverse  curvature  by  Clutter 
and  Smith  led  to  errors  in  wall  shear  stress  which  in  turn  led  to  errors 
in  transverse-curvature-induced  friction  drag  of  about  30  per  cent. 

2.  The  first-  and  second-order  treatments  of  vorticity  interaction 
were  in  substantial  agreement  over  the  forward  half  of  the  body,  and  the 
first-order  treatment  predicted  higher  wall  shear  and  thus  higher  total 
vorticity-induced  friction  drag  than  the  second-order  results  gave. 
Including  vorticity  interaction  in  the  evaluation  of  displacement  lead 

to  unrealistic  displacement  thickness  distributions  over  the  body. 

3.  The  first-  and  second-order  treatments  of  slip  and  temperature 
jump  Ted  to  significant  differences  in  slip  velocity  and  temperature  but 
small  differences  in  skin-friction  coefficient.  The  differences  in  the 
friction-drag  predictions  for  the  two  treatments  were  significant  for 
this  separate  effect,  but  the  effect  on  the  linear  combination  of  effects 
was  not  large.  A  strong  effect  of  first-order  slip  and  temperature  jump 
was  observed  when  coupled  with  all  other  higher-order  effects  at  very  low 
Reynolds  number.  However,  in  the  range  of  applicability  of  the  theory, 
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the  effects  were  small. 

4.  The  first-  and  second-order  treatments  predicted  substantially 
the  same  displacement-induced  pressure  drag;  however,  the  second-order 
treatment  predicted  substantially  higher  displacement-induced  friction 
drag.  The  difference  was  attributed  to  the  effects  of  normal  pressure 
gradient  and  external  vorticity  on  the  second-order  outer  edge  velocity 
and  temperature  since  these  effects  were  neglected  in  the  first-order 
treatment. 

5.  The  second-order  displacement-pressure  treatment  used  was  not 
unique  numerically,  but  the  error  was  small  (about  5  per  cent)  and  gave 
physically  realistic  results  for  comparison  of  this  separate  effect  with 
the  first-order  results.  Because  of  the  nonuniqueness  of  the  results, 

i  * 

only  the  sum  of  all  second-order  effects  should  be  considered  and  the 
combined  vorticity-displacement  interaction  based  upon  the  displacement 
speed  treatment  should  be  used. 

6.  Both  first-  and  second-order  treatments  substantially  over¬ 
predicted  the  experimentally  measured  pressure  distribution  for  a  sphere- 
cone  at  Mw  »  18.  These  differences  are  not  understood  since  the  compar¬ 
isons  were  made  under  conditions  where  the  second-order  theory  should  be 

applicable.  More  extensive  experimental  data  are  needed  to  clearly 

\ 

establish  the  differences,  and  further  improvements  in  the  theoretical 
model  are  needed  to  improve  the  agreement  with  experimental  observations. 

7.  The  agreement  between  first-  and  second-order  predictions  of 
Stanton  number  and  experimental  data  was  excellent  for  Mw  =  18.  Caution 
must  be  used  in  interpreting  these  results  since  limited  experimental  data 
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were  available  for  the  conditions  treated. 

8.  For  =  18  and  e  <  0.15,  good  agreement  was  found  between 
second-order  predictions  of  zero-lift  drag  and  experimental  results.  For 

e  >  0.2  at  both  a  9  and  18,  the  second-order  treatment  substantially 

.  *  .  •  ■, 

overpredicted  the  total  drag,  and  higher-order  effects  are  needed  to 

improve  agreement  between  theory  and  experiment. 

9.  For  the  conditions  treated,  the  range  of  applicability  of 
first-  and  second-order  boundary- layer  theories  is  predicted  theoret- 

t 

ically  by  the  e  corresponding  to  <S/ysh0Ck  %  1  and  from  comparison  with 
experimental  zero-lift  drag  data  was  found  to  be  e  %  0.1S. 

In  summary,  it  has  been  shown  that,  except  for  the  pressure  data 
in  the  expected  range  of  applicability  of  the  second-order  theory 
(e  <<  1),  both  first-  and  second-order  treatments  of  higher-order  boundary 
layer  effects  were  in  substantial  agreement  with  the  experimental  data. 
Beyond  the  expected  ranges  of  applicability,  the  results  from  the 
numerical  treatments  substantially  overpredicted  the  experimentally 
observed  data.  Since  boundary- layer  treatments  higher  than  second  order 
would  introduce  viscosity  into  the  outer  flow  field  and  thus  substan¬ 
tially  complicate  the  outer  as  well  as  the  inner  flow-field  solutions, 
it  is  clear  that  boundary- layer  treatments  higher  than  second  order  are 
not  the  correct  approach.  A  theoretical  model  is  needed  which  properly 
takes  into  account  all  higher-order  effects  such  as  a  viscous  "external" 
flow  and  the  effects  of  transport  properties  on  the  shock  wave.  The 
fully  viscous  shock-layer  treatment  described  by  Davis  [31]  should  pro¬ 
vide  a  valuable  tool  for  investigating  this  interesting  and  difficult 
problem. 
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APPENDIX  A 

NUMERICAL  METHODS  FOR  SOLVING  FiRST-ORDER 
BOUNDARY-LAYER  EQUATIONS 

The  treatment  of  the  first-order  boundary-layer  effects  was  based 
on  the  method  of  ciutter  and  Smith  [8]  since  considerable  work  had  pre¬ 
viously  been  done  by  Lewis  and  Whitfield  [7]  using  the  original  method 
without  vorticity  interaction  or  slip  and  temperature  jump. 

The  essential  features  of  the  method  will  first  be  desribed,  and 
then  the  modifications  required  to  include  STJ  and  vorticity  will  be 
discussed. 

The  nonsimilar  terms  in  the  momentum  and  energy  equation  were 
replaced  With  two-point  difference  expressions,  and  the  resulting  equations 
were  ordinary  differential  equations.  The  split  boundary-value  problem 
was  then  reduced  to  an  initial-value  problem  for  numerical  solution. 

•  The  procedure  used  for  solving  the  boundary- layer  equations  is 
shown  in  Figure  42.  The  momentum  equation  is  coupled  to  the  energy 
equation  through  the  fluid  properties,  and  when  the  momentum  equation 
solution  is  known,  the  energy  equation  is  a  linear  differential  equation 
in  ip  -  g  -  1.  The  method  then  was  as  follows:  First,  fluid  properties 
were  assumed  either  crudely  to  initiate  the  solution  or  from  the  solution 
at  the  previous  station  or  x  location.  Then,  the  momentum  equation  was 
integrated  from  the  wall  n  *  0  to  the  outer  edge  of  the  boundary  layer 
n  =  n^.  (The  technique  used  for  integrating  the  momentum  equation  is 
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First 

Solution 

of 

Momentum 

Equation 


Second 

Solution 

of 

Momentum 

Equation 


Fig.  42  Clutter  and  Smith's  Procedure  tor  Solving  the  Coupled  Momentum 
and  Energy  Equations  (from  Clutter  and  Smith  (8]) 
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described  below.)  Using  the  assumed  fluid  properties  and  the  resulting 
trial  solution  of  the  momentum  equation,  two  linearly  independent  solutions 
were  obtained  for  the  energy  equation.  Since  the  energy  equation  is 
linear,  a  linear  combination  of  these  two  solutions  was  obtained  satisfy* 
ing  the  wall  and  outer  edge  boundary  conditions.  New  fluid  properties 
were  obtained,  and  a  new  solution  was  obtained  for  the  energy  equation 
using  the  same  solution  of  the  momentum  equation.  This  process  was 
repeated  a  prescribed  number  of  times  or  until  two  successive  solutions 
of  the  energy  equation  agreed  within  a  prescribed  test,  that  is,  until 
<  ee  where  j  is  the  iteration  index.  Then  a  new 

w  w 

solution  of  the  momentum  equation  was  obtained  based  on  the  updated  fluid 
properties.  The  procedure  for  solving  the  system  was  repeated  a  pre¬ 
scribed  number  of  times  or  until  |  <  e^,  Typical  values 

used  in  the  calculations  without  vorticity  or  STJ  were  as  follows: 


Equation 

Iterations 

Convergence  Test 

Momentum 

5 

0.001 

Energy 

3 

0.001 

The  values  listed  above  were  input  data  and  could  be  changed  as 
required;  however,  the  machine  program  is  large  and  the  machine  time 
required  was  affected  by  the  choice  of  data.  Therefore,  the  tests  were 
made  as  loose  as  consistent  with  the  problems  solved. 

The  choice  of  method  used  to  solve  the  momentum  equation  was 
important  since  the  equation  is  nonlinear  and  under  certain  conditions 
the  solutions  were  found  not  to  converge.  Figure  43  shows  typical  trial 
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Fig.  43  Transformed  Velocity  Profiles  Resulting  from  Trial  Solutions  of 
the  First-Order  Momentum  Equation  in  Clutter  and  Smith  Variables 
(from  Clutter  and  Smith  [8]) 
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solutions  for  the  momentum  equation.  The  technique  used  is  called 
"shooting"  or  "hunting"  since  solutions  are  sought  which  satisfy  the  wall 
and  outer  boundary  conditions,  namely  ♦’CO)  ■  -1  and  4»  *  Cn^)  a  0  neglect¬ 
ing  vorticity  and  STJ. 

Briefly  the  method  used  was  as  follows:  A  trial  guess  of  the 
solution  <|>"w  was  made  (either  input  or  from  a  previously  converged 
solution).  Hie  momentum  equation  was  integrated  from  the  wall  n  ■  0  to 
n  ■  n„  or  until  |$*|  >  K  where  K  was  an  input  constant  which  was 
typically  K  »  0,5  (in  Figure  43,  K  =  1.0).  The  value  of  <|>"w  was  increased 
or  decreased  if  <J’,(nmax)  <  K  or  ^'Cn^x)  >  K,  respectively,  until  ■ 

nw.  The  procedure  was  continued  until  three  solutions  were  found,  one  of 
which  was  a  high  solution  with  K  >  >  0,  and  one  a  low  solution 

-K  <  4)'^  <  0,  and  the  third  solution  was  either  high  or  low.  The  three 
solutions  were  then  interpolated  for  a  solution  which  satisfied  the  outer 
boundary  condition  at  n  . 

The  method  used  to  integrate  both  the  momentum  and  energy  equations 
was  described  by  Clutter  and  Smith  as  a  four-point  extrapolation  and 
interpolation  formula,  and  they  gave  the  details  of  the  methods  used. 

Because  of  the  iterative  technique  used  by  Clutter  and  Smith  for 
solving  the  momentum  and  energy  equation,  it  was  easy  to  modify  the 
boundary  conditions  to  include  vorticity  and/or  STJ.  With  vorticity 
the  outer  boundary  condition  becomes 

0» 

<l>,„  *  &(*)  /  Pe/P 
o 
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and  with  STJ  the  wall  boundary  condition  becomes 

<*'w  +  1  ■  £'w  *»  [py/(p*y*)](p*/p)Fs  fn(n»0) 


r  , ' 


and 


+  1  "  gW 


,  15  py 

\  *  rix 


p*  3h  ue2  -,2] 

—  F»nr*-r  fiJ 


where  Fs  is  defined  in  Equation  86c, 

The  remaining  significant  change  occurred  in  the  density  profile 
when  vorticity  was  included.  Since  the  pressure  is  constant  across  the 
boundary  layer,  the  density  profile  becomes  (for  a  perfect  gas) 

0e  h  2(H/U„2)  -  f'2 

“  ■  2tvuei)  .  1 

Since  with  vorticity  f 1  («)  >  1,  the  boundary  layer  cannot  extend  to 
infinity  because  at  some  distance  h  ■+  0  or  T  -*•  0  and  p  Setting  the 

static  enthalpy  h  to  zero  in  the  above  equation  and  solving  for  f* 
gives  the  maximum  or  limiting  edge  velocity  possible 

f'mx  ■  l/d  '  V"e>l/2 

The  edge  of  the  boundary  layer  (nw)  must  lie  below  this  line  of  infinite 
density.  Thus  for  large  vorticity  effects  (large  ft),  the  value  of  n» 
was  reduced  as  the  solution  progressed  along  the  body.  It  was  often 
necessary  to  reduce  the  an  step  size  in  the  integration  of  the 
momentum  and  energy  equations  since  values  of  nm  <  1  were  not  uncommon. 
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In  order  to  reduce  the  step  size  An,  it  was  necessary  to  stop  the  cal¬ 
culation,  interpolate  all  quantities  which  were  functions  of  n,  and 
continue  the  calculation  with  the  new  step  size.  The  procedure  was 
straightforward  but  time  consuming,  especially  when  it  was  necessary  to 
interpolate  the  resuits  several  times  fo i  a  given  set  of  conditions. 
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APPENDIX  B 

FIRST-  AND  SECOND-ORDER  BOUNDARY-LAYER  DISPLACEMENT  EFFECTS 

In  both  first-  and  second-order  treatments  it  was  necessary  to 
obtain  a  solution  for  the  perturbed  inviscid  outer  flow.  The  outer  flow 
solution  was  obtained  over  the  geometric  body  perturbed  by  the  first- 

order  displacement  thickness  defined  as  follows: 

00 

6*  ■  /  [l  -  pu/(peue)]dy 
o 

Inverse  blunt  body  and  characteristics  solutions  were  obtained  over  the 
effective  body.  The  method  used  was  a  modification  of  one  developed  by 
Lomax  and  Inouye  [32] .  The  modification  was  necessary  to  properly  des¬ 
cribe  the  effective  body  since,  in  general,  no  simple  means  was 
available  to  describe  the  entire  body.  An  accurate  approximate  scheme 
was  developed  which  permitted  fitting  the  body  by  three  separate  curve- 
fit  procedures  with  proper  matching  conditions. 

Approximate  treatments  have  been  used  in  the  past  to  perturb  the 
outer  flow  field.  Davis  and  Fliigge-Lotz  perturbed  the  first-order 
inviscid  flow  solutions  by  assuming  the  effective  body  to  be  a  body 
approximated  by  shifting  and  expanding  the  original  geometric  body. 

This  treatment  had  the  desirable  feature  of  being  independent  of  the 
expansion  parameter  e,  and  thus  the  results  were  quite  general.  The 
disadvantage  which  ruled  out  this  method  was  that  the  shifted  and 
expanded  body  technique  was  limited  to  the  nose  region  (see  [13]). 
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Since  inviscid  solutions  were  required  over  the  entire  effective  body  and 
especially  over  the  conical  afterbody,  a  new  and  more  general  method  of 
solution  was  developed. 

A  crude  attempt  was  previously  developed  by  Lewis  and  Whitfield  [7] 
to  treat  the  displacement  effect.  In  that  treatment,  the  spherical  nose 
was  expanded  or  contracted  by  the  displacement  thickness  at  the  initial 
data  line  just  downstream  of  the  sonic  line,  and  an  approximate  curve  was 
faired  between  the  spherical  nose  and  a  curve  fit  of  the  downstream 
effective  body  data.  That  method  worked  reasonably  well  except  at  very 
low  Reynolds  numbers  in  which  case  the  displacement  thickness  became  large 
and  the  region  covered  by  the  approximate  curve  fit  became  large  in  extent. 

The  method  used  herein  was  an  improvement  on  the  scheme  of  Lewis 
and  Whitfield  [7]  in  that  the  numerical  accuracy  was  substantially 
improved.  The  perturbed  spherical  nose  was  approximated  with  small  error 
by  a  conic  where  the  bluntness  ratio  was  computed  by  fitting  the  effec¬ 
tive  body  at  the  initial  data  surface.  This  surface  was  determined  by  a 
rotationally  symmetric  line  normal  to  the  body  where  M  =  1.05  at  the 
inviscid  wall.  A  curve  was  then  fitted  to  the  effective  body  data  such 
that  r,  dr/dz,  and  d^r/dz^  were  continuous  on  the  surface  at  the  initial 
line  (M  =  1.05)  and  also  at  the  point  where  the  curve  joined  the  curve 
which  fitted  the  data  over  the  conical  afterbody.  The  intermediate 
curve  fit  was  also  required  to  have  no  inflection  points,  and  the  best 
fit  of  the  data  in  the  least-squares  sense  which  satisfied  the  other 
criteria  was  chosen. 

The  results  of  the  cuTve-fit  procedure  are  shown' in  Figure  44. 
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The  fit  became  poorer  as  values  of  e  increased  from  0.1  to  0.5  since  an 
inflection  in  the  effective  body  curve  occurred  at  the  higher  values. 
However!  in  a  least-squares  sense  the  fit  remained  good  even  though  the 
extent  of  region  increased  substantially.  Allowing  inflection  points  to 
occur  in  the  region  near  z  *  0.8  would  have  induced  a  secondary  shock  and 
substantially  altered  the  flow  field  over  the  conical  afterbody.  Experi¬ 
mental  data  did  not  indicate  any  such  secondary  shock ,  and  such  a  shock 
was  not  allowed. 

The  perturbed  inviscid  outer  flow  solution  method  (blunt  body  and 
characteristics  solution)  was  identical  for  both  first-  and  second-order 

t 

boundary- layer  treatments.  The  difference  in  application  were:  (1)  the 
way  in  which  the  outer  boundary  conditions  are  obtained.  (2)  the  effects 
of  the  normal  pressure  gradients,  and  (3)  the  effects  of  inviscid 
vorticity. 


I .  F IRST-ORDER  TR  EATMENT 

The  perturbed  inviscid  outer  flow  pressure  was  imposed  on  the  first- 
order  boundary- layer  equations,  and  consistent  with  the  first-order  theory, 
the  normal  pressure  gradient  was  neglected.  The  edge  velocity  was  obtained 
from  Bernoulli's  equation,  and  the  derivative  was  obtained  numerically. 

The  effects  of  inviscid  vorticity  were  only  implicitly  involved  through 
the  blunt  body  and  characteristics  solution.  Through  the  positive  viscous- 
induced  pressure  increment,  the  first-order  displacement  treatment  always 
produced  an  increase  in  both  pressure  and  friction  dTag.  This  result  is 


to  be  contrasted  with  the  second-order  treatment  described  below. 
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II.  SECOHD-ORDER  DISPLACEMENT  PRESSURE  TREATMENT 

The  skin-friction  coefficient  Cf  is  defined  as  follows: 

Cf„  *  [2/p;U*2]  (v*  3u*/3n*)n_Q 
■  2 e  ♦  2e^  T2(s,0) 

where 

TiCs.O)  -  (pu1n)Nb0 

t2£s,o)  -  (yu2N  ♦  y  t2  u1N)Na0 

Therefore,  if  t2(s,0)  is  independent  of  e,  then  the  friction-drag 
component  will  be  independent  of  e,  and  the  total  second-order  friction 
drag  Cd£2  a  /  Cf^  cos  a  2ir  rwds/Aj,  will  be  quadratically  dependent  on 
e»  The  first-order  friction  drag  is  given  by  C|jfl  «  e.  The  dependence  ' 
was  easily  seen  by  inspection  of  the  second-order  equations  except  in  the 
case  of  the  displacement  effect*  For  all  second-order  effects  except 
displacement,  e  did  not  appear  in  the  calculation.  However,  the  outer 
edge  boundary  conditions  for  the  second-order  displacement  pressure 
effect  depended  on  e  through  the  extrapolation  of  the  inviscid  outer 
flow  pressure  and  velocity  from  the  effective  body  to  the  geometric 
wall.  Therefore,  a  perturbed  inviscid  outer  flow  solution  was  obtained 
for  different  values  of  e  over  the  range  of  interest,  whereas  the 
other  effects  were  simply  scaled  by  from  a  single  calculation  for 
given  geometry  and  free-stream  conditions  (sec  also  Marchand,  Lewis, 
and  Davis  [18]). 

The  discontinuity  in  surface  curvature,  k  ,  caused  difficulties 
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in  numerically  evaluating  the  x-derivative  of  the  extrapolated  inviscid 
outer  flow  pressure.  The  discontinuity  was  replaced  by  an  exponential 
fairing  between  k  a  l  and  0.  The  constant  in  the  exponential  term  was 
varied  as  the  x-step  size  was  changed  until  a  smooth  transition  in  shear 
stress  was  obtained  in  the  region  of  the  sphere-cone  tangent  point.  This 
scheme  reduced  some  large  oscillations  induced  by  the  discontinuity 
without  substantially  affecting  the  total  displacement-induced  friction 
drag. 


III.  COMPARISON  OF  FIRST-  AND  SECOND-ORDER  TREATMENTS 

The  significant  difference  between  the  first-  and  second-order 
treatments  was  in  the  normal  pressure  gradient.  In  the  second-order 
theory,  the  second-order  normal  momentum  equation  was  used  to  extrapolate 
the  pressure  along  the  effective  wall  to  the  geometric  wall  or  body 
surface.  Consistent  with  the  original  first-order  boundary- layer 
equations,  the  normal  component  of  momentum  was  not  considered  in  the 
first-order  treatment,  and  the  pressure  was  assumed  constant  throughout 
the  boundary  layer  normal  to  the  surface.  Therefore,  in  the  first-order 
treatment  the  perturbed  inviscid  outer  flow  pressure  was  impressed  upon 
the  first-order  boundary  layer  as  the  edge  pressure. 

The  comparative  effects  on  the  first-  and  second-order  solutions 
were  as  follows: 

1.  Over  the  nose  where  the  longitudinal  curvature  k  *  1,  the 
pressure  on  the  surface  from  the  second-order  theory  was  either 
increased  or  decreased  from  the  effective  wall  value,  depending  on 
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whether  the  displacement  thickness  was  negative  (highly  cooled  wall)  or 
positive.  Over  the  conical  afterbody,  where  <  »  0,  the  first-  and 
second-order  treatments  yielded  the  same  surface  pressure,  assuming 
equivalent  displacement  thicknesses.  However,  because  of  the  effects  of 
the  inviscid  vorticity  along  the  wall,  the  velocity  at  the  surface  was 
always  less  over  the  conical  afterbody  based  upon  the  second-order 
theory.  This  effect  due  to  inviscid  vorticity  was  independent  of 
inclusion  or  exclusion  of  the  vorticity  effect  in  the  second-order  boundary- 
layer  treatment.  Again,  the  coupling  between  displacement  and  vorticity 
is  seen,  but  this  time  through  the  inviscid  flow  fields, 

I 

2.  The  first-order  displacement  treatment  always  induced  an 
increase  in  friction  and  pressure  drags.  The  second-order  effects  could 
be  either  positive  or  negative,  depending  on  the  sign  and  magnitude  of 
displacement  thickness,  and  the  effects  of  longitudinal  curvature  and 
inviscid  external  vorticity. 

These  two  major  factors  can  cause  the  first-  and  second-order 
displacement  treatments  to  yield  substantially  different  effects  on  such 
global  quantities  as  drag  and  total  heat  transfer;  however,  for  the 
length  of  the  cone  and  other  conditions  treated  here,  the  first-  and 
second-order  treatments  predicted  induced-drag  effects  in  good  agreement. 
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APPENDIX  C 

EFFECTS  OF  VORTICITYON  DISPLACEMENT  THICKNESS 
Assuming  the  following  form  for  the  inviscid  velocity  profile, 

ui  ■  ue  c1  *  ui  y) 

where  Uj  =  (f1  (•)  -  l)/y£«0*  then  the  definition  of  displacement  thick¬ 
ness 

5*  oo  m 

j  oju^r^dy  »  /  p^u^r^dy  -  J  purJdy 

0  0  0 

leads  to  the  following  integral  equation: 

V  1  i '' 

6v 

B  /  (Pi/Pgjtl  +  U!  y)(rw  +  y  .cos  ct)J  dy  . 
o 

oo 

=  B  /  (Pi/PgHl  +  U'  y)(rw  +  y  cos  a)3  dy 
o 

oo 

-  /  f '  Crw  +  y  cos  a)J  dn 
o 

where  B  =  [Re*  (ue/uj/x] l^2  pe/P*. 

When  the  vorticity  is  neglected,  in  *  0  and  the  integral  equation  for 
6*  reduces  to  the  classical  axisymmetric  form  of  the  equation  with  TVC, 
Results  are  shown  in  Chapter  IV  for  the  variation  of  6*  over  the  body. 
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